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ABSTRACT 

We evaluate the prompt observational signatures of the merger between a massive close- 
in planet (a 'hot Jupiter') and its host star, events with an estimated Galactic rate of ~ 0.1 - 1 
yr . Depending on the ratio of the mean density of the planet p p to that of the star p*, a 
planet-star merger results in three possible outcomes. If p p /p* > 5, then the planet directly 
plunges below the stellar atmosphere before being disrupted by tidal forces. The dissipation of 
orbital energy creates a hot wake behind the planet, producing a EUV/soft X-ray transient that 
increases in brightness and temperature as the planet sinks below the stellar surface. The peak 
luminosity Leuv/x < 10 36 erg s _1 is achieved weeks to months prior to merger, after which 
the stellar surface is enshrouded by an outflow driven by the merger. The final stages of the 
inspiral are accompanied by an optical transient powered by the recombination of hydrogen 
in the outflow, which peaks at a luminosity ~ 10 37 - 10 38 erg s _1 on a timescale ~days. 

If the star is instead significantly denser (p p /p* < 5), then the planet overflows its Roche 
Lobe above the stellar surface. For p p /p* < 1 mass transfer is stable, resulting in the planet 
being accreted on the relatively slow timescale set by tidal dissipation. However, for an in- 
termediate density range 1 < p p /p* < 5 mass transfer may instead be unstable, resulting 
in the dynamical disruption of the planet into an accretion disk around the star. Outflows 
from the super-Eddington accretion disk power an optical transient with a peak luminosity 
~ 10 37 - 10 38 erg s _1 and characteristic duration ~week-months. Emission from the disk 
itself becomes visible once the accretion rate decreases below the Eddington rate, resulting in 
a bolometric brightening and shift of the spectral peak to UV wavelengths. Optical transients 
from both direct-impact merger and tidal-disruption events in some ways resemble classical 
novae, but can be distinguished by their higher ejecta mass and lower velocity -hundreds km 
s _1 , and by hard pre- and post-cursor emission, respectively. The most promising search strat- 
egy is with combined surveys of nearby massive galaxies (e.g. M31) at optical, UV, and X-ray 
wavelengths with cadences from days to months. 
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1 INTRODUCTION 

Theoretical models of planet formation suggest that gas giants form 
beyon d the ice line at radii > several AU from their central star 
(e.g. Ilnaba et all [2003). It is thus a mystery how some exoplan- 
ets, colloquially known as 'hot Jupiters', are transported to their 
current locations at radii < 0.5 AU (e.g. Ilda & Lin! |2004|) . The 
misalignment between the orbital planes of some hot Jupiters and 
the spin axes of their host stars (Tri aud et al. 2010; Schlaufman 
l20ld - |Winn et alj|2010h suggests that not all migra te via interac- 
tion w ith the proto-stellar disk alone (although see iFoucart & Lail 
1201 lh . At least some hot Jupiters appear to have mig rated via other 
mechanisms, such as planet-planet scattering (e.g. iRasio & Ford 
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Fabr vckv & Tremaine ■|2007tlSocrates et al.l201 ]]). In both cases the 
planet initially approaches small radii on an eccentric trajectory, be- 
fore tidal dissipation circularizes the orbit. 

Guill ochonetalj|20T 1), on the other hand, argue that the cur- 
rent locations of many close-in exoplanets are inconsistent with 
migration resulting solely due to the inward scattering from orig- 
inal orbits exterior to the ice line, since the small pericenter dis- 
tances implied by such trajectories would have resulted either in 
tidal-disruption of the planet or have completely ejected it from the 
system. One interpretation of their result is that some planets have 
migrated further after their orbit is circularized, due to tidal dissi- 
pation within the star. 

The angular momentum of most hot Jupiters is sufficiently low 
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that no sta te of tidal equilibrium exists (th e system is |Parwinll880l 
unstable' ; lHutll 198 ll ; iLevrard et al.ll2009l) . suggesting that the ulti- 
mate 'end state' of ti dal dissipation is a m erger between the planet 
and its host star (e.g. I Jackson et al.ll2008l) . Depending on the qual- 
ity factor of tidal dissipation within the star Q' t ~ 10 s , the semi- 
major axes of several known hot Jupiters are sufficiently small that 
a merger wi ll indeed occu r on a relatively short timescale < 10-100 
Myr (e.g. Ili et alj[2010l) . The lack of old planetary systems with 
very short orb ital periods hints tha t mergers indeed result from tidal 
orbital decay dJackson et ai]|2009l) . 

One mi ght expect that since ~ 1 per cent of stars host hot 
Jupiters (e.g. IMavor & Ouelozll2012l). and since the Galactic star- 
formation rate is ~ 1 - 10 yr ~' jNaab & Ostriker!l2006h . then the 
Galactic rate of planet-star mergers should be < 0. 1 yr" 1 . However, 
the actual merger rate of planets with main sequence stars can in 
principle be significantly higher ~ 1 yr~' , depending on the value of 
Q' t and the rate at which hot Jupit ers are replenished b y migration 
from the outer stellar system (e.g. ISocrates et al.ll201lt see §0. A 
comparable or greater number of mergers may occur in evolved 
stellar systems, since the timescale for tidal migration T p \ unge oc R~ 5 
(eq. (T|) depends sensitively on stellar radius R+, which increases 
as the star evolves off the main sequence. In fact, irrespective of 
the efficiency of tidal dissipation, it seems inevitable that almost all 
planets with main-sequence orbital separation < 1 AU eventually 
merger with their stars, since stars expand dramatically during the 
red giant and asymptotic giant branch phases of evolution (although 
post-main-sequence mergers are significantly less energetic events 
than their main sequence brethren). 

In this paper we explore the direct observational consequences 
of the merger of gas giant planets with their host stars. Depend- 
ing on the properties of the planet-star system, we show that the 
merger event can be violent, resulting a bright and long-lived elec- 
tromagnetic signature. Given the rapid pace of recent advances in 
our knowledge of exoplanets, and the advent of sensitive wide-field 
surveys across the electromagnetic spectrum, now is an optimal 
time to evaluate the transient signatures of planet-star mergers and 
how best to go about detecting them. 

This paper is organized as follows. In ^2]we summarize the 
properties of known hot Jupiters and use them to estimate the rate of 
planet-star mergers in galaxies similar to the Milky Way. In ^3]we 
identify three qualitatively different merger outcomes, which are 
distinguished in part by whether the planet fills its Roche radius at 
an orbital separation a t which is larger or smaller, respectively, than 
the stellar radius R+ . If a, < R+ then the planet is not disrupted until 
it has fallen below the stellar photosphere. In Sj4]we describe the dy- 
namics of the planetary inspiral and argue that such 'direct-impact' 
merger events are accompanied by a extreme ultraviolet (EUV)/soft 
X-ray transient that originates from the hot gas created behind the 
planet as it penetrates the stellar atmosphere Q4.4. U . This emis- 
sion, which ceases weeks-months prior to merger, is immediately 
followed by rising optical emission powered by the recombination 
of hydrogen in an outflow from the stellar surface, which peaks on 
a timescale ~days after the planet plunges below the stellar surface 

If a, > R*, then the planet instead overflows its Roche radius 
above the stellar surface. If mass transfer is stable, then the planet 
is consumed on the relatively long timescale set by the rate of tidal 
dissipation; this results in a relatively short-lived mass-transferring 
binary, but is unlikely to produce a bright electromagnetic dis- 
play. If, on the other hand, mass transfer is unstable (or if the ini- 
tial planet trajectory is eccentric), then the planet is disrupted into 
an accretion disk just outside the stellar surface ('tidal-disruption 



event'). Thermal radiation from the disk, and from super-Ed dingt on 
outflows from its surface, powers a bright optical-UV transient 
(iJS). In ijfj] we discuss our results, including a concise summary 
of the predicted transients from planet-star mergers Q6.lt ; detec- 
tion prospects and strategies, in particular with surveys of nearby 
galaxies at optical, UV, and X-ray wavelengths ( 36. 2b ; and a brief 
conclusion, which focuses on unresolved theoretical issues and fu- 
ture work ( 36. 3t . 

Our work here builds on, and is broadly consistent with, past 
theoretical work into the observational s ignatures of binary stel - 
lar mergers (e.g. iTvlenda & S oker 20061 ISoker & Tvlendall2006l) . 
eve nts whic h Soker & Tvlenda ( 2006) refer to as 'mergebursts ' (see 
also lBearetalj201lh . 



2 GALACTIC RATE OF PLANET-STAR MERGERS 

Dissipation of tides raised on a planet and on its host star can 
act to change their orbital separation. Two important categories of 
tidal interactions can be termed "circularization tides" and "syn- 
chronization tides." Dissipation of the former kind acts to damp 
eccentricity at constant orbital angular momentum £.. Since £. oc 
Va(l - e 2 ), where a is semi-major axis and e is eccentricity, cir- 
cularization of the orbit (e — » 0) at constant £. shrinks the binary 
separation. 

Dissipation of synchronization tides instead occurs due to 
asynchronous rotation of the bodies with respect to the orbital mean 
motion, even in the absence of eccentricity. In particular, if a planet 
orbits faster than its host star rotates, then tidal torques transfer an- 
gula r momentum from the planet to the star, le ading to orbital de- 
cay dDarwiJl88dlHutlll98ll ; ISoker & Tvlendall2006l) . The charac- 
teristic timescale for this process (the tidal 'plunging time' T plungc ) 
is found by integratin g the tidal dissipation equations assu ming a 
slowly rotating star dLevrard et al.l2009t,lBrown et al.bOlTb . which 
can be estimated as (see also Appendix IB1 

Tp, u „ge ~ n ^ (1. 

where Q\ = 3Q*/k2 is the modified tidal quality factor of the star 
(<2* is the specific tidal dissipation function and k 2 is the Love num- 
ber); n = 2n/P ^ (GMt,/a 3 ) 112 is the orbital mean motion; P is 
the orbital period; a is the orbital semi-major axis; R+ is the stellar 
radius; and and M p are the stellar and planetary masses, respec- 
tively, where the latter is normalized to the mass of Jupiter Mj. 

Figure [TJ shows the cumulative number of planets as a func- 
tion of tidal plunge time, calculated by applying equation (l) 
to the distribution of ~ 160 known transiting extrasolar planets 
dSchneideretailbOUh FI This Tpi un ge distribution can be used to 
estimate the rate of planet-star mergers if one can extrapolate the 
known sample of exoplanets to those of the galaxy as a whole. In 
order to make this conversion, we assume that the observational 
selection function is given simply by the geometric probability of 
transit (R*/a). We also assume that a fraction ~ 10~ 5 of all stars 
in the Milky Way have been searched by transit surveys (G. Bakos, 
private communication), such that each transiting planet used in 
Figure[T]"counts" as 10 5 (a/i?*) planets in the Galaxy. 

Dashed lines in Figure Q] represent the expected plunge- time 

1 http://exoplanet.eu 
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distribution if one assumes that mergers occur with a steady-state 
rate of 10 6 /2* Y r 1 an d 0.5(10 6 /<2'*) yr respectively. A compari- 
son of the observed distribution to these models indicate that, at low 
values of plunge time ~ 10 6_7 ((2i/10 6 ) yrs, roughly one merger 
occurs in the Galaxy every 2((2' lk /10 6 ) years. Obvious deviations of 
the cumulative rate from the model predictions occur at large values 
Tpiunge- These indicate the onset of additional effects not included in 
our simple model, such as a dep endence of Q' t on orbital period 
(e.g. lGoodman & Dicksot]|[l998l) ; the onset of additional selection 
effects; or a break-down of the steady-state assumption. The latter 
could result from the o ngoing injection of planets from larger radii 
(e.g. lSocrates et al.l201 1 ), or by insitu pl anet formation following a 
binary star merger (e.g. Martin et al.ll201 ll) . Although our estimate 
of the merger rate from Figure [TJis uncertain by at least an order 
of magnitude, the fiducial range of values ~ 0.1 — 1 yr" 1 that we 
find are still sufficiently high that transient surveys with duration 
comparable to this mean interval ~year — decade could in princi- 
ple detect such events within the Milky Way or nearby galaxies, 
provided that the resulting emission is sufficiently bright. 

As mentioned in the Introduction, the rate of planet-star merg- 
ers could be enhanced as a result of post main sequence evolution. 
After the star evolves off of the main sequence, its larger radius R t 
significantly decreases the plunging timescale Tpi^ge oc R t 5 (eq. 171 : 
even in the subgiant phase, for which R* ~ few R Q ), more than com- 
pensating for the shorter duration of post main sequence evolution. 
This suggests that a sizable fraction of all stars with semi-major 
axes within a < 0. 1 AU may be consumed JSato et al.l20 08; Hansen 
|2010| ; iLlovduOlU) . Since the fraction of stars w ith planets in this 
radius range approaches unity dMavor et al.l201 IB . the Galactic rate 
of planet-star mergers may approach a significant fraction of the 
star-formation rate ~ 1 - 10 yr" 1 . 

Mergers are likely to be even more common during the 
gian t phase when the st ellar radius expands to R+ ~AU 
(e.g. lSiess & Lividll999al lbl; lRetter & Maromll2003l) . The low den- 
sity of the stellar envelope in this case would almost certainly re- 
sult in what we define as a direct-impact merger (Sj4]l. However, the 
much lower gravitational energy released near the stellar surface 
in such an event as compared to the merger with a dwarf or sub- 
giant results in much fainter transient emission that is more chal- 
lenging to detect ( although these events may have other observ- 
able consequences; Retter & Maroml 12003; Nordhaus et al.ll201ol ; 
Nord haus et alj|2oTlh . For this reason we focus on planet mergers 
with compact stars, even though their total rate may be somewhat 
lower. 



3 THREE POSSIBLE MERGER OUTCOMES 

So long as no mass transfer occurs between the planet and the star, 
tidal forces control the rate at which the orbit decays. As the planet 
inspirals towards the star, its fate depends on whether it fills its 
Roche-lobe above or below the stellar surface. In the latter case the 
planet plunges directly into the stellar atmosphere. This produces 
an event that we define as a 'direct-impact' merger event, the re- 
sult of which we describe in detail in Sj4] If Roche-lobe overflow 
instead begins above the stellar surface, then the planet is either 
dynamically disrupted into a compact disk around the star ('tidal- 
disruption' merger event; ^5}, or it gradually transfers its mass over 
a much longer timescale (stable mass transfer). 

The volume of the Roche-lobe of the planet in a synchronous 
orbit about a star of much greater mass is approximately given by 
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Figure 1. Cumulative number of planet-star mergers in the Galaxy (solid 
blue line) as a function of plunging time r p i ungi . (eq. (TJ), the latter normal- 
ized to a fiducial value Q' t ~ 10 6 for the rate of tidal dissipation in the 
star. Overlaid dashed lines indicate the expected distributions for steady- 
state merger rates of (10 6 /Q' t ) yr~' (red) and 0.5(10 6 /Qi) yr~' (green), 
respectively. Deviations of the model predictions from the observed distri- 
bution at larger values of T p i ungc indicate the presence of additional effects 
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Figure 2. Sample of 160 transiting extra-solar planets in the space of mean 
stellar density p* versus mean planetary density p p . The plunge time of 
each planet Tpi^g,, (eq. (TJ) is denoted by the size and fill of each marker. 
The shape of each marker indicates the predicted fate of each system as 
described in and Appendix [A] assuming that the mass and radii of the 
planet and star are the same upon merger as their current values. Black 
triangles indicate planets that directly impact the stellar surface (Pp/p* > 5; 
'direct-impact' merger events; 5|4); green squares to planets that undergo 
stable mass transfer above the surface (Pp/p* < 1); and red circles to planets 
that undergo unstable mass transfer and are shredded into an accretion disk 
around the star (1 < p p /p* < 5; 'tidal-disruption events'; 5f5). 



V RL - 0.5(M p /MJfl 3 toai&Blandfor3l201lh . Roche-lobe over- 
flow thus begins when the orbital semi-major axis reaches the value 



a, - 2fl*(p,/p p ) 1/3 , 



(2) 



where p\ = 3M*/4jtRI and p p = 3M p /4nR 3 p are the mean density 
of the star and planet, respectively. 

At the orbital separation given by equation ((2) the Ll La- 
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Table 1. Outcomes of planet-star mergers based on known transiting planets 



Planet Sample 


Merger Outcome 


Fraction (%) 


Tphmge < (<2;/10 6 )10 7 yrs 


Direct-Impact 


33 


- 


Tidal-Disruption 


22 




O LdUlfJ IVlclj 3 1 l i. i 1 I.M L. 1 


44 


T p i U nge<(G'J10 6 )10 lu yrs 


Direct-Impact 


13 




Tidal-Disruption 


25 




Stable Mass Transfer 


62 


All Systems 


Direct-Impact 


12 




Tidal-Disruption 


22 




Stable Mass Transfer 


66 



grange point is located at a radial distance X, ^ 0.7(M p /M lr ) 1/3 a, ^ 
lARp from the center of mass of the planet. The condition that 
Roche-lobe overflow occurs entirely below the stellar surface is 
thus given by a, < R+ + X,, or 

2(pjp p ) >13 < 1 + 1.4(tf p /tf*). (3) 

If equation 10 is satisfied, then a direct-impact merger occurs (ij4j. 

If Roche-lobe overflow instead occurs outside the star (i.e. if 
eq. [3] is not satisfied), then the evolution of the system is more 
subtle. In this case mass transfer is stable(unstable) when mass loss 
from the planet results in the Roche-lobe growing faster(slower) 
than the radius of the planet. Here we describe the qualitative dis- 
tinction between the system evolution in these two cases, the details 
of which are provided in Appendix lAl 

If overflow occurs when the orbital separation is sufficiently 
large, then material crossing LI has sufficient specific angular mo- 
mentum to form a disk well above the stellar surface. Viscous re- 
distribution of angular momentum within the disk causes most of 
its mass to accrete onto the star, whereas the majority of the angu- 
lar momentum is transported to large radii where it can be trans- 
ferred back into the orbit via tides exerted by the planet. Since 
mass is transferred from the less massive planet to the more mas- 
sive star, this (approximate) conservation of orbital angular mo- 
mentum acts to widen the orbital separation. Since the Roche-lobe 
volume V RL oc a 3 increases rapidly as the orbit expands, mass trans- 
fer is temporarily halted. Orbital decay due to tidal dissipation then 
again drives the system together, resulting in a series of overflow 
episodes that slowly exhaust the planetary material on the tidal de- 
cay timescale T p i lmei: (a = a, > R+) > 10 3 yr (eq. Q]). This is the 
stable mass transfer case. 

If, on the other hand, Roche-lobe overflow occurs only at a 
sufficiently small planet-star separation, then material crossing LI 
plunges directly onto, or circularizes just above, the stellar surface. 
In this case the bulk of the mass and angular momentum of the 
accreting material is deposited into the stellar envelope, which is 
capable of accepting it due to its considerable inertia. Since or- 
bital angular momentum is lost to the star, the semi-major axis of 
the orbit (and hence the Roche volume of the planet) remains un- 
changed, whereas the planet expands upon mass loss (Appendix 
lAl . This results in unstable mass transfer and hence to dynamical 
tidal-disruption of the planet on a characteristic timescale of several 
orbital periods ~hours (^5). 

Figure|2]shows the sample of transiting planets used in Figure 
[JJin the space of mean stellar density p* versus mean planetary 
density p p . Each planet is marked according to both the value of its 
tidal plunging time (eq. IT)) and by our best estimate of its ultimate 



fate upon merger ('direct-impact'; 'stable mass transfer'; or 'tidal- 
disruption'). The criteria for the latter are evaluated assuming that 
the current masses and radii of the planets and their host stars are 
identical to those at the time of merger (the validity of which we 
discuss in j 36, 3t . We find that if the density ratio is sufficiently high 
(Pp/p* >5), then the the planet is disrupted below the stellar surface 
(direct-impact mergers). Planets somewhat denser than their host 
star (1 < Pp/p* < 5) will instead be dynamically disrupted close to 
(but outside) the stellar surface (tidal-disruption). Finally, the least 
dense planets (p p /p* < 1) typically undergo stable mass transfer. 

Table [TJ summarizes the percentage of known transiting plan- 
etary systems that are predicted to experience various outcomes 
upon merger. Out of the entire sample of planets in Figure|2] we es- 
timate that ~ 12 per cent undergo a direct-impact merger. This frac- 
tion, however, rises to ~ 33 per cent (3 out of 9) when considering 
just those systems with the shortest lifetimes r p i un g e < (<2'*/10 6 )10 7 
yrs. We estimate that ~ 1/4 of systems will undergo unstable mass 
transfer and tidal-disruption, while the remaining (approximately 
half to two thirds) are consumed less violently via stable mass 
transfer. Note again that our estimates are implicitly restricted to 
the merger of planets with main sequence stars. If one also includes 
mergers during the post main sequence, then the expansion of the 
star decreases the value of p* and hence would increase the fraction 
of planets that undergo direct-impact or tidal-disruption. 



4 DIRECT-IMPACT MERGER 

If the planet is disrupted well inside the star (p p > 5p*), then most 
of the gravitational energy released during this process, £ 6rav > 
GM+Mp/Rt = 4 x 10 46 (M p /Mj) ergs, goes into heating the stel- 
lar interior. This additional energy will most likely act t o enhance 
the stellar luminosity (however, see |Podsiadlowskill2003l) . but only 
on the cooling timescale of the envelope at the disruption radius. In 
the Sun, for example, the radiative cooling timescale at r ~ 0.5R B 
(where the mean interior density is p ~ 10 g ctrT 3 ) is^| t coo i ~ 10 6 
yr. Disruption of a Jupiter mass planet would thus increase the stel- 
lar luminosity by an amount L ~ E irm /t coa \ > 10 33 erg s _1 , com- 
parable to the normal solar luminosity. Thus one characteristic of 
direct-impact mergers is a long-term change in the stella r evolu- 
tion and appearance dPodsiadlow ski 2003; Tvlenda 2005). Unfor- 
tunately, the amplitude of this variation is typically modest and its 
duration is typically longer than is accessible to transient surveys. 

In this section we instead describe the more prompt signature 
that originates as the planet first interacts with the external layers 
of the star, before disappearing into the interior. The planet impacts 
the stellar surface on a grazing orbit with the Keplerian velocity 
\\ = (GM+ZR+y 12 . This leaves a hot trail of shocked stellar gas be- 
hind it with a characteristic virial temperature T vlI ~ 10 s — 10 7 
K, resulting in a powerful UV/X-ray transient for a duration of 
weeks-months (we show below that the spectrum peaks at a some- 
what lower temperature than r v i r ). In what follows we assume that 
the orbit of the planet is nearly circular at the time of merger, which, 
as we show in Appendix [B] is justified for the sample of known 
transiting extrasolar planets. However, in ^6] we briefly describe 
how our conclusions change if the planet instead impacts the star 
on a highly eccentric orbit. 

2 This may be an overestimate of the true cooling timescale, since heating 
from the disruption will increase the entropy near the point of disruption, 
thereby altering the stellar structure. 
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Figure 3. Schematic diagram of the direct-impact of a planet with the stellar surface (5j4). Show in top and side views is the geometry of the hot layer left 
behind as the planet grazes the stellar atmosphere (yellow). The latitudinal width of the wake perpendicular to the planet l± is set by the width of the planet 
at the depth of the stellar scale height H. The length of the wake parallel to the planetary motion is set by the distance l\\ = \\t zao \ that the planet passes in a 
cooling time f coo i, with a typical distance of several scale heights. An (uncertain) fraction e lac i < 1 of the thermal power E generated as the planet sinks through 
the atmosphere (eq. 0) is radiated directly from the surface (typically EUV/soft X-rays; see eq. (9)), while the rest of the energy (likely the majority) goes 
into driving an outflow further downstream (see Fig.|5]for a schematic diagram of the outflow geometry). 



4.1 Inspiral dynamics 

In this section we describe the dynamics of the inspiral of the planet 
through the stellar atmosphere. Then below we address the proper- 
ties of the heated stellar surface ( 94.2b and the resulting electro- 
magnetic emission from the merger ( 94.4b . 

As the planet first encounters the stellar atmosphere of density 
p*, it loses angular momentum and energy via turbulent drag, caus- 
ing the planet to slowly spiral inwards with time. The drag force per 
unit area is given by / d = Cap* vjJ/2, where C<i is the dimensionless 
drag coefficient, which is of order unity for a supersonic flow with a 
high Reynolds number. The atmospheric density decreases rapidly 
with radius, with a (unperturbed) scale height H ~ 10~ 3 - lO~ 2 R t 
in the outer layers which is much less than the planetary radius 
R p ~ OAR*. The torque on the planet is thus dominated by the 
fluid that intercepts the planet within a distance « H from its inner 
edge closest to the star. Simple geometr^ shows that in the limit 
that H <k R p , the effective cross section of the planet for angular 
momentum loss is approximately given by A p ss H x L ± , where 
L x =! R p 2 H 112 is the characteristic width of the star perpendicular 
to the direction of motion at the penetration depth H (see Fig. [3}. 
This expression is accurate for H < R p , as is valid through the final 
stages of the inspiral. 

In calculating the torque on the planet, we are justified in us- 
ing an unperturbed value for the atmospheric scale height because 
the cooling timescale of the hot gas behind the planet is typically 
shorter than the orbital period (see below), such that the stellar at- 
mosphere should have approximately returned to its original state 
by the time that the planet has passed around the star. We are also 
justified in neglecting the effects of the gravitational field of the 

3 We neglect gravitational focusing, since this only increases the effective 
cross section of the planet by a factor ~ M p R+/M+R p <K 1. 



planet on the stellar scale height since the gravity of the star typi- 
cally dominates near the stellar surface. Even if the gravity of the 
planet is important, the planet-crossing timescale ~ R p /\\ ~ 10 2 s 
is generally shorter than the timescale required for the atmosphere 
to adjust thermally ~ H/c s ~ 10 3 (/f/ 10 _2 i? o ) s. Likewise, we ne- 
glect the effects o f nonaxisymmetric gravita tional tides raised by 
the passing planet dRicker & Taaml2008ll2012l) . assuming that they 
are negligible compared to other forms of drag. 

The drag torque faA p x R+ extracts angular momentum from 
the planet at the rate / = M p v^v r /2, resulting in an inward velocity 
given by 

Vr * c^rI'V^m; 1 ^- (4) 

One may define the characteristic 'infall' time t- m llw from any ra- 
dius as that required for the planet to migrate inwards by a scale 
height, since this is the distance to encounter denser gas (and hence 
substantially stronger drag): 

where f orb = 2nR+/\\ ss 10 4 s is the orbital time at the stellar sur- 
face and in the second expression (and hereafter) we assume that 
R+ = R , R p = 0AR o , and Cd = 1. Equation l[5} shows that the 
infall time is much longer than an orbital time in the outer lay- 
ers of the star where the stellar density is low, but that the rate of 
infall increases with time as the planet plunges deeper, such that 
f in,drag < forb (dynamical plunge) for p* > 0.01 - 0.1 g cm 4 . Infall 
is clearly an exponential runaway process. 

The stellar density above which inward migration due to gas 
drag dominates over that due to tidal dissipation can be estimated 
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by comparing ?j n _d mg to the corresponding infall time due to tidal 
migration 



H 



<Wud -U T#»vL** ~ 17 yr . 10 _ 2 ^ q M 1Q6 / 1 m, 



(6) 



where the factor (H/R+) corrects for the fact that T p i unge is defined 
in equation (TJ as the migration time across a distance of order the 
orbital separation a x R+.By equating $5$ and ((6} using a realistic 
profile H(p+) corresponding to a main sequence solar type star (see 
below), we find that gas drag dominates over tidal drag (assuming 
<2* = 10 6 ) at stellar densities greater than p* > 10 -5 (10~ 35 ) g cm -3 
in the case of a planet of mass M p = l(10)M o . In both cases this 
first occurs at a characteristic infall time of fj n _d rag ~ fin.tidai ~ 100 
days. 

The rate at which energy is dissipated by the inward migration 
of the planet is given by 

GM,M p v,. 



2Rl 



* 6xl0 37 ergs~ 



2x lO^ergs" 1 



J gcm 



H 



\Q- 2 Rr 



3/2 



H 



lQ- 2 R n / \ M 



I orb 



(7) 



where we have used equation ((4) for v, . In the last expression we 
have substituted the infall time t m ljmg from equation ((5}, so this ex- 
pression is valid only when f^drag < 'in.tidai- 

Figure [4] shows several quantities relevant to the dynamics of 
the inspiral of a hot Jupiter planet into a solar mass star, as a func- 
tion of the increasing stellar density p t encountered as the planet 
moves inwards. These include the density scale height of the stellar 
atmosphere H\ the dissipated power E (eq. (7J); and the planet infall 
time t m = min[f ln ,tidai,fin,dm g ] (eqs. (UU) assuming = 10 6 , nor- 
malized to the orbital time r or b ~ 10 4 s. For t m two cases are shown, 
corresponding to different planet masses M p = lMj and 10Mj, re- 
spectively. The times at which the inward migration rate due to gas 
drag exceeds that due to tidal dissipation (i.e. t iri 

,drag ^ ^in.tidal) 

marked in the top panel with vertical dotted red lines. Our calcu- 
lations employ the stellar structure H(p+) of a 1M Q ZAMS star at 
a characteristic age t = 5 G yr, as calculated using the EZ stellar 
evolution code (Paxton 2004). 

Figure [4] shows that as the stellar density increases from ~ 
l(T 5 (10- 4 )g citT 3 to ~ lO^lO" 1 ) g cnT 3 , the infall time t m de- 
creases from ~ 10 3 / or b ~100 days to dynamical infall ~ f orb ~hrs, 
for planetary masses M p = IMj(IOMj), respectively. The dissi- 
pated power increases over this timescale from ~ 10 33 erg s -1 
to ~ 10 40 erg s -1 in the M p = lMj case. The evolution in the 
M p = 10Mj case is similar, except that the final power at the point 
of dynamical infall reaches a higher value ~ 10 41 erg s due to the 
greater gravitational energy released (eq. Q). Although the rate of 
dissipated power becomes very high during the final stages of the 
merger, below we show that E represents only an upper limit on the 
radiated electromagnetic power. 



4.2 Properties of the Hot Post-Planetary Layer 

As the planet moves through the stellar atmosphere, the dissipation 
of shocks and turbulence creates a wake of hot gas behind it (see 
Fig.[3]for a schematic illustration). The temperature of this hot gas 
r h is determined by equating its specific thermal energy ~ c 2 with 
the power generated E/M ~ (Cd/2)vJ; per rate of swept-up stellar 
mass M = App+xy, where c s is the sound speed in the hot gas. The 



evolution of is shown in the bottom panel of Figure [4] indicates 
that 7h increases from ~ 10 6 K to ~ 10 7 K with increasing p*, 
i.e. much hotter than the surface layers of the unperturbed star T ~ 
10 4 - 10 5 K. 

In order to calculate the emission from the stellar surface we 
must first estimate the size of the hot post-planetary layer. The 
width of the layer perpendicular to the planetary motion l ± ss 
(RpH) >/2 is set by the width of the planet at the depth of the at- 
mospheric scale height. The length of the layer along the direction 
of motion Zy ~ f coo iv k is instead set by the rate at which matter cools 
downstream, where f coo i is the cooling timescale. 

At early times, when the planet encounters the low-density 
outer layers of the star, the cooling timescale is set by the rate 
at which radiation diffuses vertically / d iif ~ 3(///c)t, where t = 
K m p+H is the vertical optical depth and K es is the opacity, the dom- 
inant form of which is Thomson scattering given the high temper- 
atures T h > 10 6 K of the shocked gas. Although heating may sub- 
stantially increase the scale height of the atmosphere further down- 
stream from the planet (see below), we are justified in adopting an 
unperturbed value for the stellar scale height H in calculating f diff 
provided that cooling by diffusion is much faster than the rate of 
thermal expansion of the heated material (i.e. if t&g <K H/c s ), as is 
indeed satisfied when r is sufficiently small. 

As the depth of the planet (and hence r and t$g) increases, 
convection becomes more efficient than radiation at transporting 
energy generated by the planet to the stellar surface. If convection 
dominates the cooling, then the length of the hot layer Zy can in- 
stead be determined by equating the dissipated power E ss 
with the maximum luminosity that can be carried by convection 
£<c,max ~ (ph v cmax/2) x C±'||)> where Ph is the density of the post- 
planetary gas and v craax = a c c s (where a c < 1) is the maximum 
convective velocity. For the moment we assume that ph ~ p*, 
i.e. that the gas in the hot layer does not have time to expand ap- 
preciably from its unperturbed value (and hence is drastically out of 
pressure equilibrium; we discuss the validity of this assumption be- 
low). In this case the length of the hot layer is Zy ~ H/a\ ~ few H, 
corresponding to a convective cooling timescale f c00 i = f c>vert ~few 
xH/vt. 

Combining results, the emitting area of the hot post-planetary 
layer is approximately given by 



A„ a l ± X Zy « Rl l2 H ll2 v k t cooh 



(8) 



where f coo i ss min[rdiff = 3//r cs /c, t CrVett = 3///vk]. The transition 
between radiative and convective cooling (tm > r c ,vcn) occurs at 
stellar densities p* ~ 1(T 5 g cirT 3 and is marked by the vertical 
dotted line in the bottom panel of Figure|4] 

Our estimate of f c . vcrt above assumes that hot gas behind the 
planet does not expand appreciably from its unperturbed value, re- 
sulting in a hot wake of length Zy ~few xH when f d iff ^ ? c,ven- 
In reality, the full extent of the heated region may be significantly 
larger than this estimate because the cooling timescale ? c vcrt oc p^ 1 
increases as the atmosphere expands, a process which itself occurs 
on a timescale ~ H/c s similar to the cooling time. Furthermore, in 
34.3| below we show that the hot post-planetary gas powers an out- 
flow with a velocity ~ c s ~ Vk, which itself acts to cool the gas on 
a similar timescale ~ H/c s . 

Because the timescales for convective cooling / C-V eit, atmo- 



4 The cooling and radiative diffusion timescales are equal because the en- 
ergy density in the post-planetary wake is dominated by radiation. 
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10 _ 6 1Q _5 1Q _4 10 - 3 1Q -2 1Q -1 

Stellar Density p„ (g cm -3 ) 



Figure 4. Quantities relevant to the dynamics of the inspiral of a hot Jupiter planet into a solar mass star in the case of a direct-impact merger ($4), as a 
function of the increasing stellar density p* encountered as the planet moves inwards. Quantities shown include the stellar density scale height H (dotted blue 
line, bottom panel); the dissipated power E (eq. (7); dark solid line, top panel); and the planet infall time fj n , in units of the orbital time f or |, « 10 4 s (eqs. |5|6| ; 
dot-dashed line, top panel). For ?;„ two cases are shown, corresponding to different planet masses M p = lMj (purple) and lOMj (orange), respectively. Times 
after which the inward migration timescale due to gas drag tj aj j at (eq. (5)) is less than that due to tidal dissipation fj n ,tidal (eq- 0, assuming Q' t = 10 6 ) are 
shown in the top panel with vertical dotted red lines. Also shown on the bottom panel are the temperature Th of the hot layer on the stellar surface behind the 
passing planet (dark solid line) and the effective temperature T e ff of radiation from the hot layer (orange dot-dashed line). The stellar structure H(p ir ) used in 
our calculation is that of a 1M Q ZAMS star of age t = 5 Gyr (see text for details). 
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spheric expansion, and wind cooling are of the same order|f| it 
seems reasonable to assume that a modest fraction e ra j of the gener- 
ated thermal energy will transported to the stellar surface and radi- 
ated just a few scale heights behind the planet (see Fig.O, resulting 
a luminosity Leu\/x = tmtiE. Below, we adopt a fiducial value of 
e md = 0.1, since obtaining a more accurate estimate of e rad will 
require a detailed radiation-hydrodynamical simulation of the in- 
teraction between the planet and stellar atmosphere, a task beyond 
the scope of this paper. In what follows below it should be kept in 
mind that the value of e ra[ ] is highly uncertain (by at least an order 
of magnitude). 

The temperature of the emission from the hot post-planetary 
layer is approximately given by 

1/4 



\A/,cr sb 



5 x lCf K 



10- 3 gcirr 3 



H 



OMR, 



-1/8 



(9) 



where cr s b is the Stefan-Boltzmann constant and we have used 
equation ([8) for A n . The evolution of r eff is shown in the bot- 
tom panel of Figure [4] assuming e rad = 0.1. If the spectrum is ap- 
proximately that of a blackbody, then the peak photon energy is 
£peak ~ 3A;r cff = 260(r cff /10 6 K) eV. Note that T cfs increases from 
an initial value of ~few 10 5 K (extreme ultraviolet; EUV), to a value 
T eff ~ 10 6 K (soft X-rays) at later times closer to merger. In 34.41 



' eff 

below we describe under what conditions this radiation escapes the 
vicinity of the star and hence is actually observable. 




V„,~ V v 



4.3 Outflows 



Figure 5. Schematic diagram of the outflow driven from the stellar surface 
by the energy generated as the planet grazes the stellar atmosphere. The 
mass loss rate of the wind M w is comparable to the rate at which mass is 
swept up by the planet (eq. 1111 ). whereas the asymptotic velocity of the 
wind is assumed to be similar to the stellar escape speed ~ \\. On small ra- 
dial scales the wind is non- stationary since it is launched from a (time vary- 
ing) region well-localized behind the passing planet. However, on larger 
scales r > R m ~ v w totb ~ 6Rq (where r orD is the planet orbital period) 
the outflow is reasonably well approximated as a spherically symmetric, 
steady-state wind. 



An outflow is likely from the hot post-planetary layer because the 
total generated power ~ E (eq. (7)) greatly exceeds the (local) Ed- 
dington luminosity L cd d « GM^A^c/R^k^, with their ratio given 
by 



E 

LvAA 



2A h ~ 



0.5/f e! pA(vk/c)*2x 10 4 | 



10~ 3 g cm - 



(10) 



where we have used equations © and ([8} for v r and Ah, respec- 
tively. Equation l llOt shows that E/L c m > 10 2 - 10 5 for the relevant 
range of stellar densities p* ~ 10~ 5 - 10~ 2 g cirT 3 encountered 
during inspiral. 

If the fraction 1 - e ra( j ~ 0(1) of the thermal power that is 
not radiated immediately behind the planet is used to power a wind 
further downstream, then the resulting energy per unit mass > c 2 
is similar to that required to escape the potential well of the star 
A fraction of the thermal energy will go into exp anding a 



hot envelope (e.g. lTvlendall2005l : ISoker & Tvlendall2006l) , while a 
comparable fraction will go into the kinetic energy of an outflow. 
Assuming that the asymptotic velocity of the wind v w is comparable 
to the escape speed of the star ~ Vk, then the mass loss rate M w 



In principle cooling may also occur non-radiatively due to the redistri- 
bution of heat latitudinally by convective motions or subsonic excavation 
flows (e.g. Melosh 1989), since in this case thermal energy is lost to PdV 
work via sideways expansion. However, the condition for latitudinal trans- 
port to dominate over vertical transport is that l± > H, i.e. R p > H, which 
is in practice only satisfied when the atmospheric density is sufficiently 
high pi, > 10~ 2 g cirT 3 that the planet is already sinking dynamically 

(i.e. £i n ,drag ~ ^orb)- 



of the outflow will be comparable to the rate at which the stellar 
atmosphere is swept up by the planet, i.e. 

\3/2 



A p p*v k a 7x 10- 2 gs 



10- 3 gcm- 3 /\l0- 2 «, 



H 



(11) 



Note that M w is many orders of magnitude greater than the normal 
mass loss rate in the solar wind ~ 10 12 g s _1 . 

Since the bulk of the wind is launched from a narrow layer 
behind the passing planet, the outflow will be non- stationary and 
highly asymmetric near the stellar surface. However, on larger 
scales, greater than the distance /fj„ « v w ? OI t ~ 6R that wind ma- 
terial traverses in an orbital time, the flow will develop a structure 
similar to that of a spherically symmetric, steady-state wind. Fig- 
ure|5]shows a schematic illustration of the wind geometry, showing 
an asymmetric inner cavity (r <k R m ) which connects onto a steady 
state wind at larger radii r » R m . At sufficiently large radii, the 
density profile of the wind is approximately given by 



47T?' 2 V„ 



7x 10- 



H 



io- 2 ^ 



3/2 



r 

6R, 



(r > Rm), 



(12) 



where we have used equation Jilt for M w . 

If the material in the wind expands adiabatically, then its tem- 
perature T v (r) will be a factor (p w /p*y smaller than the surface 
temperature of the hot layer r h , where /? ~ 1/3 - 2/3 for a gas 
with adiabatic index T ~ 4/3 - 5/3 (depending on whether the out- 
flow is dominated by gas or radiation pressure). Near the inner edge 
r = ~ 6/f Q of the steady-state wind, typical temperatures are in 
the range T w ss 10 3 5 - 10 s 5 K forp* ~ 10~ 5 - 10~ 2 g cirT 3 . 
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Figure 6. Properties of the radiation from a direct-impact planet-star merger, as a function of time until the dynamical plunge ("merger") of the planet below 
the stellar surface fj n- drag (e<3- 0)- Top and bottom panels show cases corresponding to planet masses M p = lMj and lOMj, respectively. Quantities shown 
include the total rate at which energy is dissipated by the sinking planet E in units of 10 30 erg s (solid black line; eq. (7)); EUV/X-ray luminosity from the 
stellar surface just behind the hot planet in units of 10 30 erg s -1 Leuv/x = £ ra d£ (dashed grey line; see 34.4. U assuming e ra d = 0.1; peak photon energy from 
the hot layer behind the planet on the stellar surface £„eak (orange line; in eV); optical depth through the wind surrounding the star to electron scattering T es 
(dashed purple line) and to free-free absorption T a t, s- ff (at energies E a £ pca k and accounting for the additional photon path-length due to electron scattering; 
dashed green line). Vertical red dashed lines mark the times after which the photon luminosity from the stellar surface Leuv/x is suppressed due to inelastic 
electron scattering (r cs > 40; T a b s ,cs > 1) and free-free absorption (T a b s ,ff > 1). A dashed blue line shows an estimate of the optical luminosity from the 
photosphere of the optically-thick stellar wind in units of 10 30 erg s ( 34.4.21 . which is estimated as L opl = min[L rec , L pea k], where Lrec is the luminosity 
due to hydrogen recombination in the outflow (eq. 1151 ) and L pca k (eq. 1171 ) is the peak luminosity achieved following the final inspiral stage, after which the 
outflowing material better resembles a dynamically ejected shell than a steady-state wind. 
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4.4 Electromagnetic Radiation 

4.4.1 EUV/X-ray Emission from the Early Inspired 

Figure [6] shows the power dissipated during inspiral for planet 
masses M p =1 and 10 Mj cases (top and bottom panels, respec- 
tively), similar to that shown in Figure |4]but now plotted as a func- 
tion of time until merger t m (eq. (5j), starting at the time / ~ 100 
days after which the infall rate due to gaseous drag exceeds that 
due to tidal dissipation. Shown with a dark gray line is the (orbit- 
averaged) luminosity L E uv/x = e miE from the hot post-planetary 
layer, assuming e la d = 0.1 (see above). Shown in the bottom panel 
with an orange line is the peak photon energy L pea k = 3£L eff 
(eq. (9)) in eV. 

Radiation from the stellar surface is in general not isotropic 
due to the geometry of the planet with respect to the hot radiat- 
ing layer (see Fig. [3}- Rather, for a typical observer, emission is 
beamed into a fraction = Af2 rad /47r ~ 0.25(H/R p ) [/2 < 0.1 of 
the sky, as determined by the solid angle A£l rai j ~ n(H/R p ) ,/2 be- 
tween the stellar surface and the planet. Provided that radiation is 
free to propagate from the stellar surface to infinity unattenuated, 
one predicts that the observed emission will vary periodically on 
a timescale t olb ~ hours, with an emission duty cycle f b ~ 10 per 
cent and a peak luminosity L E uv/x/./6 ~ 10L E uv/x- The shape of the 
light curve will depend both on the viewing angle of the observer 
and on the precise geometry of the radiating layer. 

If radiation indeed escapes from the stellar surface unat- 
tenuated, one would predict that Leuv/x should increase from 
~ 10 34 - 10 36 erg s -1 at t ~100 days prior to merger, up to 
L E uv/x ~ 10 4 ' - 10 42 erg s~' in the final hours before merger. In re- 
ality, however, such a high luminosity during the final stages of the 
inspiral is not actually achieved because outflows from the stellar 
surface ( 34.3b block the escaping radiation once the mass loss rate 
becomes sufficiently high. In particular, even if radiation escapes 
the vicinity of the stellar surface immediately behind the passing 
planet unhindered, it must necessarily pass through the mass lost 
from previous orbits (Figure|5). Since this "orbit- averaged" outflow 
is first encountered on a radial scale > R m ~ 6R Q , the characteristic 
optical depth through the wind is r ~ p w R m K w , where p w (eq. 1121 ) 
and k w are the density and opacity of the wind, respectively, at radii 
r ~ R m . 

The opacity k w depends on the ionization state of the wind 
a radii r > Rj n , As described above, the mean energy of radia- 
tion from the stellar surface E^.^ > 100 eV is sufficiently high 
> 13.6 eV to ionize hydrogen. Whether the wind is fully ionized 
thus depends on the ratio of the emission rate of ionizing pho- 



Nft, « fieOftV,), x n e r , viz. 
No, 



JVn. vo.i/ 



ll0 4 K/ 



H 



-3/2 



r 

6R, 



10" 3 gcm- 3 / \10 4 K/ \10- 2 i? Q 

where Ry = 13.6 eV is the Rydberg; « c ss p w /m p and vj, = 
(fcL w /m e ) 1/2 are the number density and thermal velocity of elec- 
trons in the wind, respectively; and cr^, =s 2 x 10~ 21 (L W /10 4 K) 
cm 2 is the radiative recombination cross section (Milne relation), 
where L w is normalized to a typical value (see discussion after 
eq. 1121 ). Equation J 1 3b shows that Nmn/Nfb is » 1 for all rele- 
vant values of the stellar density and wind temperatures, indicating 
that the wind will be fully ionized at all radii r > R in as long as out- 
flow is transparent to ionizing radiation, the conditions for which 
we now describe. 

The two dominant sources of opacity in fully ionized mat- 
ter are electron scattering and free-free absoiption. The electron 



scattering opacity k cs ss 0.4 cm 2 g~' is independent of frequency, 
whereas the free-free opacity k s cc p w 7\!/ 2 v~ 3 on radial scale « R in 
is given by 



K ff u„, - io-wg->m ^r 3/i 



\Tpeak 
/4 



10- 3 scm- 3 / V 10 



m 

h \" /8 / r w v 
- 2 rJ \10 4 K/ 



(14) 



where we have normalized the incident frequency v to the peak 
frequency v pi:ak = E^.^/h ss 3kT s . fl /h of the emission from the stel- 
lar surface using equation (|9). The low characteristic value of Kg 
in equation J14t shows that, at the peak frequency, electron scat- 
tering dominates the opacity by several orders of magnitude (al- 
though free-free opacity becomes dominant at lower frequencies in 
the Rayleigh- Jeans tail). 

Figure|6]also shows the optical depth through the wind to elec- 
tron scattering (dashed purple line). When t cs < 1, correspond- 
ing to t m > 100 days, radiation escapes freely from the stellar 
surface, resulting in a periodic light curve with a peak luminosity 
> 10Leuv/x as described above. On the other hand, at later times 
when t cs » 1 scattering acts to effectively "smear out" this light 
curve variability. Nevertheless, at lest initially the observed lumi- 
nosity will continue to equal the orbit-averaged value Leuv/x even 
when t cs > 1 , because electron scattering is (to first order) energy 
conserving. 

As the time until merger decreases, the wind becomes increas- 
ingly opaque, such that radiation from the stellar surface is com- 
pletely absorbed. The two dominant processes which thermalize 
photons in the wind are (1) inelasticity due to repeated electron 
scattering; and (2) free-free absorption. Since the number of scat- 
tering events experienced by a photon is ~ r 2 s (for t cs > 1), while 
a fraction f sc ~ hv w ^lm t c 2 ~ 5 x 10~ 4 of the original energy is 
lost per scattering, inelastic electron scattering results in apprecia- 



ble energy loss once t cs > f sc 



1/2 



~ 40. Figures [6] shows 



that for M p = l(10)Mj this condition (marked by a vertical dashed 
line) is satisfied at times t m < 7(70) days, when the EUV/X-ray 
luminosity is L EU v/x ~ (e ra d/0.1)10 36 erg s" 1 . 

Surface radiation is also attenuated by free-free absorption, 
which becomes more effective than suggested by equation d 14b due 
to the longer path length created by multiple electron scatterings. 
One can account for this by defining an effect ive mean free path for 
thermal absorption as T abs ff = T es (*- ff //f es ) 1/2 tevbicki & Lightmanl 
1 19791) . where Kg is the free-free opacity at the peak frequency 
(eq. 1141 ). A second vertical dashed line Figure [6] marks the in- 
fall time below which free-free absorption is important (T abs ff > 1), 
indicating that this process is somewhat less effective than inelas- 
(13) tic electron scattering at thermalizing radiation near the peak of the 
EUV/X-ray spectrum. 

High-energy photons absorbed by the plasma are re-emitted as 
lower energy photons, for which the the free-free opacity is higher 
(e.g. Kg oc y~ 3 ; eq. 1141 ). Once the timescale for photons to dif- 
fuse out of the wind becomes longer than the expansion timescale 
~ r/v w , most of this energy is lost to PdV work. Directly observ- 
able emission from the stellar surface is thus be substantially less 
than Leuv/x = fradL after times for which r abs = min(T abs , cs ,r abs , ff ) 
> 1, as Figure [6] shows is satisfied in the weeks-months prior to 
merger. The EUV/X-ray luminosity should thus peak at a value 
LEuv/x( T abs ~ 1) ~ ( e iad/0.1)10 36 erg s -1 , after which the emission 
is strongly suppressed as the region surrounding the star becomes 
increasingly enshrouded by the merger-driven outflow. 
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4.4.2 Optical Emission from the Final Inspiral 

In the previous section we showed that hard emission from the 
stellar surface is suppressed by absorption long before the final 
merger. However, lower frequency emission from the outflow it- 
self becomes visible, and continues to strengthen, as the outflow 
power approaches its peak value E ~ 10 40 - 10 4 ' erg s _I . Most of 
the initial thermal energy in the wind at the stellar surface is lost to 
PdV expansion before it can be radiated at the photosphere of the 
wind. However, the recombinatio n of hydrogen contributes an a ddi- 
tional source of luminosity (e.g. lKasen & Raimrez-RuizhoiOl) and 
acts to slow the rate at which the temperature decreases (effectively 
reducing the adiabatic index). The energy released by hydrogen re- 
combination, if efficiently radiated, produces a 'luminosity' given 
by 



(M„/ra p ) x Ry 



Ry . 



0.03£. 



(15) 



Recombination is possible in the outflow only once energetic pho- 
tons from the stellar surface are no longer able to keep the wind 
ionized; Figure [6] shows that this condition (r abs > 1) occurs 
on a timescale of weeks-months prior to merger, when the out- 
flow power E ~ 10 37 erg s~' and the recombination luminosity is 
L rcc ~ 3 x 10 35 ergs" 1 . 

In ^5] we calculate the detailed properties of photospheric 
emission from outflows with very similar properties to those con- 
sidered here (j j4.3t . but powered instead by super-Eddington accre- 
tion following the tidal-disruption of a planet (see Figs.[8]|9]l. These 
calculations show that the radiated luminosity is indeed a signifi- 
cant fraction of the recombination luminosity ~ L lcc (eq. 1151 ). and 
that the radiation peaks at optical wavelengths (photospheric tem- 
perature ~ 4 - 6 x 10 3 K). 

Our results in §5\ show that equation l !15t provides a reason- 
able estimate of the optical luminosity during the earliest stages of 
the final inspiral t iulimg ~ weeks-months. However, equation l |15t 
overestimates the luminosity when the infall approaches its dynam- 
ical stage and E becomes very high. This is because the assumption 
that the outflow is in steady state is no longer valid once the plane- 
tary infall time ?in,d ra g (eq. (5)) becomes shorter than the expansion 
timescale of the wind f exp = R mc /vi, at the recombination radius 
R Ttc . Since this condition f cxp > fi n> d rag is satisfied on a timescale ~ 
day-week prior to merger, equation (115) is only a valid estimate 
of the optical luminosity until the optical luminosity L opl reaches a 
value- 10 37 - 10 38 ergs" 1 . 

A significant fraction of the energy released, and total mass 
ejected, during the merger occurs during the final few orbits before 
the dynamical plunge. An accurate treatment of the emission dur- 



ing this final phase when t a 



requires modeling the ejected 



mass as that of a discrete 'shell' of ionized material with a mas^| 
M c j « M p /10. The luminosity under this type of evolution does 
not peak on the same (short) timescale ~ f or b of the final plunge 
(i.e. when E peaks), but rather on the longer timescale set by when 
the expanding ejecta becomes sufficiently cool for hydrogen to re- 
combine and radiate its energy. 

We may crudely estimate the peak luminosity L pca k and peak 
timescale ? pcak of this final dynamical phase by equating f peak ~ 
£rec/£peak with the timescale ? cxp ~ i? P hoto/vk required for the ejecta 
to expand to a sufficiently large radius ^ p hoto for the energy released 



This factor 1/10 results because the planet falls a distance comparable 
to its own radius ~ R p ~ O.IR* in the final plunge, resulting in sufficient 
gravitational energy released to eject ~ 1/10 of its own mass. 



by recombination to be radiated, i.e. t pc .^L pcli ^ ss £ rcc , where L peak ~ 
47T<r s bT' 4 c .ft photo is the thermal luminosity, £ rec ~ 0.1(M p /ra p )Ry is 
the total energy released by recombination, and r rec « 6000 K is 
the temperature of the photosphere (set by the recombination of hy- 
drogen; see Fig.[8}. This results in the following analytic estimates: 

/ r rec \- 4/3 /m p \ 1/3 



J peak : 



„ ,/ r rcc \ 4/3 /m p \ 2/3 

2 x 10 37 ergs-' — — — 
S V6000K/ \Mj 



(17) 



Thus, for M p ~ 1 - 10Mj the final merger is accompanied by an 
optical transient of luminosity ~ 10 37 - 10 38 erg s -1 and character- 
istic duration ~days (extending to times after the merger itself is 
complete). As we describe further in j 36.2l this event may resemble 
a classical nova, although is distinguished by its higher ejecta mass, 
lower ejecta velocity, and shorter typical duration. In Figure [6] we 
plot the optical luminosity with a dashed blue line, assuming that 
L opt ss L rec at early times (when the approximation of a steady-state 
wind is valid), but then capping L opl ss L peaK at late times when the 
evolution of the ejecta is dynamical. 

FigurefTJshows a schematic diagram summarizing the various 
phases of light curve evolution for the case of direct-impact plane- 
tary mergers. 



5 TIDAL-DISRUPTION EVENT 

If the planet undergoes Roche Lobe overflow above the stellar 
photosphere, then mass transfer may become unstable (ij3j see 
Fig. 0. In this case the planet is completely consumed in just a 
few orbital perio ds, resulting in a tidal-disruption event (see also 
iBear et al.ll20Hl) . Since the Keplerian velocity greatly exceeds the 
sound speed of the planetary material, a shocked virialized torus is 
formed around the star with an angular momentum approximately 
equal to that of the original orbit. A fraction of the disrupted planet 
will impact and mix with the stellar surface, whereas the majority 
will remain above the surface in rotationally-support configuration 
with characteristic radius R d0 ~ a t > R Q (in some wa ys similar 
to the 'excretion disks' observed around Be stars; e.g. iLee et al.l 
fl99lh . 

The high temperature T > 10 6 K in the disk implies that it will 
be fully io nized. Turbulence gener ated by the magnetorotational 
instability dBalbus & Hawlevlll998l) . or by shear instabilities in the 
layer between the disk and star, will then cause matter to accrete 
inwards at a characteristic rate 

M ~ — p - 



5 \0.1 )\M,\M \R e 



-3/2 



HJr\ 
0.5 ] 



(18) 



where 



R\ 



R, 



v \0.1J Mj \R a 



0.5 



(19) 



is the viscous accretion timescale (typically ~days-weeks); v = 
ac s H A is the effective kinematic viscosity; c s = H^Sl is the mid- 
plane sound speed in vertical hydrostatic equilibrium; Hi is the 
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EUV / X-Ray Emission (Stellar Surface) Optical Emission 

j (Outflow Recombination) 




- months-yr - wks-months + day(s) + ~1 Myr 



Time with Respect to Merger 

Figure 7. Schematic light curve of a direct-impact planet-star merger (SgJ. At early times, EUV/X -ray emission (blue) originates from the hot layer behind 
the planet as it grazes the stellar atmosphere (see Fig. [5). Initially the signal is periodic on the orbital timescale f or |, ~hours with a duty cycle ~ 10 per cent 
set by the narrow solid angle through which radiation escapes the cavity between the planet and the stellar surface. The rate of atmospheric heating increases 
as the planet sinks deeper, creating an outflow from the stellar surface with a mass loss rate that increases with time (Fig. [5). Electron scattering wipes out 
the X-ray periodicity months to a year prior to the merger, once the optical depth through the wind r cs > 1 (dashed vertical line). Scattered (~isotropic) high 
energy radiation continues to escape for weeks-month prior to merger, before being entirely suppressed by other processes (inelastic electron scattering and 
free-free absorption) weeks-month prior to merger (r a b s > 1; dashed vertical line). At late times, thermal optical emission from the photosphere of the wind 
(red) becomes increasingly bright, due in large part to the energy released by the recombination of hydrogen. Although initially the optical luminosity traces 
the rate of mass loss from the wind, the final dynamical stages of the merger in effect results in the ejection of a shell of ionized material from the surface. 
Thermal emission from this shell powers an optical transient lasting ~day(s) after the merger, somewhat similar to a classical nova (see eqs. 1161 . 1171 ). Finally, 
gravitational energy released by the planet deep inside the star is radiated on the much longer Kelvin-Helmholz cooling timescale of the stellar envelope ~ 10 6 
yr (green). 



vertical scale height of the disk; fl = (GMj/r 3 ) 1 ' 2 is the Keplerian 
orbital frequency; and a ~ 0.01 - 0.1 is the parametrized dimen- 
sionless viscosity dShakura & Sunvaevlll973l) . 

Equation shows that the initial accretion rate M is typ- 
ically ~ 10 2 - 10 4 times higher than the Eddington accretion rate 
Medd,* - 7 x W 12 (M i ,/M e )(R i ,/R e ) g s~' at the stellar surface, de- 
pending on the values of M p and a. At such high accretion rates 
the disk cannot cool efficiently and is thus geometrically thick with 
Hi ~ Rd/2. This is because the timescale for photons to diffuse out 
through the midplane 



fd.ff - (3HJc)T d * 6 x 10 7 s ffl |^ 



(20) 



as/ v 0.5 

is much longer than the accretion timescale f v ; sc (eq. 1191 ) over 
which gravitational energy is released, where Td a piH A K cs is the 
vertical optical depth, p d = M v l2nH A R 1 i is the midplane density, 
and k cs is the electron scattering opacity. More generally, the con- 
dition tag > f v j sc is satisfied at radii less than the critical "trapping 
radius" 



R« - —. — - R* 



M 



4nc "\M Edd .*y 
For a thick disk with Hi 



(21) 



Rd that evolves due to the vis- 



cous redistribution of angular momentum, the accr etion rate onto 
the star evolves wi th time approximately as (e.g. |Pringlelll98ll; 
iMetzger et al]|2008h 



M(t). 



M , 

M (f/wr 4/3 , 



t < t v , 



(22) 



If the evolution of the disk conserves total angular momentum^] 
Ji oc MiiGM+Rif 12 , where M A and R A are the total mass and 
(mass-)averaged disk radius, then the outer edge of the disc vis- 
cously spreads outwards in time as ^d K Mi 1 , i.e. 



R*, 

RAt/K, sc ) 213 , 



t > t v „ 



(23) 



Equations i22\ and d23b do not strictly apply once the outer edge 
of the disk is able to cool efficiently into a thin disk (i.e. once 



7 In reality, a fraction of the original angular momentum of the orbit will 
go into the outer layers of the star rather than into the disk. Nevertheless, 
the self-similar evolution of the disk radius R A oc f 2 ' 3 should still obtain at 
late times, although the normalization of R A may be somewhat smaller than 
our estimate in equation (23). 



© 2012 RAS, MNRAS 000,[H?? 



Transients from Planet-Star Mergers 1 3 



10 



0.1 - 



0.01 



1 1 1 


1 '. 


\ 7cd i Pw/ 10 " 10 g cm ~' 








■ 


: 
■ 


lor. I 


* - 


■ ^" / 


ii ; 


' ' ' 





0.1 



10 



10' 



10 J 



10 4 



10- 



10 4 



10- 



10' 







: L rec (10 M erg s"') 


— -^C 


V J : 


■ T,„ (K) N 






^s. || i 

\l (10 w erg s" 1 ) h | 



10 

time (days) 



100 



Figure 8. Example solution showing the structure of an outflow powered 
by super-Eddington accretion following the tidal-disruption of a Jupiter- 
mass planet, as a function of the optical depth from the surface r w . The 
solution is shown at a time t = f v j sc a 1.5 days, when the outflow velocity 
is v w sa 400 km s and the radiated luminosity is L ss 1.2 X 10 37 erg s 
(see Figure [9). The disk is assumed to have viscosity a = 0.1 and initial 
radius Rjq = 1.5i? . Properties shown include the wind temperature 7" w 
(black line); wind density p w (red line); ionized fraction of hydrogen X, 
(blue line); and effective adiabatic index y a( j (orange line). Also shown with 
black vertical dashed lines are the wind trapping radius R u w = r(r w = 
c/3v w ) and the photospheric radius /? p h = r(r w = 2/3). 



/? d > R tI ). These expression nevertheless provide a reasonable ap- 
proximation toJhe_viscous evolution predicted by more realistic 
models (e.g. lCannizzo et alJll990l )Fl 

The timeline of the disk evolution can thus be summarized as 
follows. At early times after the disk forms, photons are trapped 
throughout the entire disk, i.e. R tI » Ri (M » Medd,*)- Since 
R tt /Ri <*■ t~ 2 , as time proceeds the outer edge of the disk evolves 
to become 'untrapped' and radiatively efficient (i.e. R A > R„) on a 
timescale ~ (Mo/M Ed d,*) l/2 f V isc ~month when R u ~ R A » 10/? v . 
After this time, the trapping radius continues to move inwards, 
reaching the stellar surface at a time 



(a) M„ = lMj 



Mo 



3/4 



80days(Ar(^f 4 (^) 7/ \ (24) 
y \0.l) \Mj) \M S J 

At this point the outer radius of the disk is located at 

-^edd = Redi(t = 4dd) 

Thus, depending on the value of a ~ 0.01 - 0.1 and M p ~ 1 — 10Mj 
we expect f c dd ~month-few years and R^ w 10 - 300R o . 

In the following sections we estimate the properties of the 
thermal radiation from the accretion disk at times t < f cd d Q5.lt 
and t > f edd (HJ2j. 



5.1 Emission from Super-Eddington Outflows (t < t sii ) 

At early times t < f e dd when the accretion rate is highly 
super-Eddington, powerful outflows are likely to occur driven 



For a geometrically thin a-disk dominated by gas pressure and Kramers 
opacity, we find that M oc t 25 and Sd K ' at times t » ? v ; sc . 
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Figure 9. Properties of thermal emission from accretion disk outflows, as 
a function of time since planetary disruption, during the super-Eddington 
accretion phase (t < £gdd; eq. 1241 ). Top and bottom panels show cases cor- 
responding to a disrupted planet of mass M p = 1 Mj and 1 OMj , respectively. 
The wind mass loss rate is assumed to be a fraction / w = 0. 1 of the accre- 
tion rate M (eq. 1221 ). while the thermal energy at the base of the wind is 
assumed to be a fraction e t d = 0.1 of the outflow kinetic energy. Quanti- 
ties show include the radiated bolometric luminosity in units of 10 33 erg 
s (solid blue line); effective temperature of radiation at the photosphere 
r e (f = r w (r w = 2/3) (solid red line); and temperature of the wind at the 
"advection" radius a T w (r w = c/3v w ) (solid red line). Shown for 
comparison is the maximum luminosity due to the recombination of hy- 
drogen L rec (eq. 1151 ). Note that at early times the radiated luminosity is a 
significant fraction of L rec , while at late times L > L rcc once the photosphere 
recedes to near the outflow surface, such that thermal energy advected from 
the base of the wind contributes most of the luminosity. 



by radiation pressure, either directly from the accretion disk 
(e.g. ldhsuga et alJ2005h or from the stellar surface, which is heated 
by energy released at the shear interface between the disk and 
star (as in the case of d irect merger; !Q. Numerical simulations 
(e.g. lOhsugaeTal]|2005h and observations of super- Eddington ac- 
creting systems (e.g. SS433: iBegelman et alj|200r3) both suggest 
that a substantial fraction of the accreting mass is lost in outflows, 
such that the wind mass loss rate is M w ss / W M with / w ~ 0(1). 
If one assumes that the total gravitational energy released as plan- 
etary material settles from the point of disruption to the stellar sur- 
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face goes into the kinetic energy of the outflow, we estimate that 
the maximum value of / w is ~ 0.3 - 0.5 for planet-to-star density 
ratios in the range p p /p+ ~ 1 — 5 corresponding to unstable mass 
transfer; hereafter we adopt a fiducial value of /„ = 0.1. 

Because the outflow is highly optically thick, emission from 
the accretion disk itself is not visible at these early times. Neverthe- 
less, bright thermal emission can originat e from the photosphere of 
the wind, as has been stud ied before (e.g. IRossi & Begelmanl2009l ; 
IStrubbe & Ouataerti r2009) in the context of other super-Eddington 
systems. In this section we calculate the light curve and tempera- 
ture of the emission from these super-Eddington outflows. 

We assume that the wind is launched from a characteristic 
radius R w (t) = min[2?tr, i?<j] ~ 1 - 10^ Q set by the maximum ra- 
dius where the disk is locally super-Eddington, where R tI (t) and 
Ra(t) are calculated using equations d21t and J23) , respectively. 
We furthermore assume that the asymptotic speed of the wind ap- 
proximately equals the escape speed from the launching radius 
v w (r) x (GM+/R t ) [/1 « 100 - 400 km s~', as is a common propert y 
of thermal pressure-driven winds (e.g. lLamers & Cassinelli|[l999h . 
At most times the accretion rate (and hence the wind mass loss 
rate) evolves on a timescale that is relatively slow compared to the 
time since disruption (see eq. 1221 ) or the expansion timescale at 
the photosphere. In this case it is reasonable to approximate the 
density profile of the outflow as that of a steady-state wind: 



Pw(r,f) : 



M w (r) 
4m- 2 v w 



exp[-(r/r cdgc ) 2 ] 



(26) 



where the exponential cut-off accounts for the outer edge of the 
ejecta at r ~ r c a K = v w t. The optical depth through the wind to 
radius r at any time t is given by 



r„(r, r) 



J p w K w dr', 



(27) 



where K w (p w , T w ) is the Rosseland mean opacity and T w (r) is the 
temperature profile of the wind. 

At small radii in the wind (high t„.) we calculate T v (r) by 
assuming that entropy is conserved from the base of the wind at 
r = R w , out to the wind 'trapping' radius R tIVI above which radi- 
ation freely escapes without further loss of energy]^ The latter is 
defined as the radius above which the expansion timescale in the 
wind / exp ss r/v w is greater than the timescale for outwards radia- 
tive diffusion fdiff.w ~ 3r w r/c. The thermal energy at the base of the 
wind is assumed to be a fraction e th < 1 of the wind kinetic energy 
p w v 2 /2. The value of e th is highly uncertain but should be of order 
unity if the outflow represents matter preferentially heated to virial 
temperatures, i.e. above the entropy of the midplane. We assume a 
value of e th = 0. 1 in our calculations below, but our results do not 
depend sensitively on its precise value. 

Outside the trapping radius the wind is effectively stationary 
on the expansion timescale and can thus be approximated as a ra- 
diative stellar atm osphere. For r > R tIM the t emperature gradient is 
thus given by (e.g. lHansen & KawaleJI 19941) 

dr 



3 L 



16 Anr 2 cr sb T 3 



(r > R KVI ). 



(28) 



Here L is the photon luminosity of the wind. Its initial value at the 
trapping radius is set by the advection of thermal energy across the 
trap ping radius, i.e. L(r = R a v/ ) = (167r/3)<r >b r*r z | r ^ w x (v w /c) 
(e.g. lRossi & Begelmanl2009T) . At larger radii L gains contributions 
from the additional energy input from hydrogen recombination, i.e. 



dL Kl . = (M w /m p )|dX,|xRy where \dXj\ is the change in the ionized 
fraction Xj{r). 

We calculate the ionization state Xj{r) of the wind as- 
suming local thermody n amic equilibrium (Saha equation), as in 
IKasen & Ramirez-Ruij i20ld) . The opacity of the ejecta k w is 
determined using data from the Opacity ProjecQ assuming the 
metallicity of the gas to be five times solar (a reasonable approxi- 
mation to the composition of gaseous giant planets). The observed 
luminosity L(t) = 4nr 2 <r s i,T 4 and effective temperature T c (f(0 = T w 
of the radiation are those near the photosphere of the outflow 
(t w = 2/3). Note that since the trapping radius is a function of 
r w (r), which depends on the temperature structure of the wind, the 
self-consistent structure of the wind must be solved by iteratively 
"guessing" the location of the trapping radius. 

Figure [9] shows our calculation of the time evolution of the 
luminosity L; outflow temperature at the trapping radius T Rliv = 
T w (r = r a d v ); and effective temperature r e g = (L/4na- s \,R 2 h )' 14 of 
the emission, where R ph = r(r w = 2/3) is the radius of the photo- 
sphere. In the top and bottom panels, respectively, we show cases 
corresponding to different assumed planetary masses M p = lMj 
and M p = 10Mj. In all calculations we assume a wind mass loss 
parameter f v =0.1; disk viscosity a = 0.1; and initial disk radius 
Ra = l.5R , resulting in an Eddington time of f c( id ~ 80(450) days 
for M p = l(10)Mj (eq. 1241 ). For comparison we also show the 
maximum luminosity due to hydrogen recombination L lcc (eq. 1151 ) 
with a dashed blue line. 

Figure [9] shows that the wind luminosity peaks at a luminos- 
ity L ~ 10 37 (10 38 ) erg s -1 for M p = l(10)Mj, in both cases on 
a timescale ~ t visc ~ 1 day, before declining with time approx- 
imately as L oc f -1 (somewhat more shallow than M oc r 4/3 ). 
During this time the effective temperature is relatively constant at 
r eff x 5000 - 7000 K, due largely to the sensitive temperature de- 
pendence of the opacity in this temperature range. At early times 
(/ < 10 days) most of the luminosity results from the energy re- 
leased by hydrogen recombination, i.e. L < L rec (eq. 1151 ), as op- 
posed to the initial thermal energy advected from the base of the 
wind. This implies that the evolution near peak luminosity is rela- 
tively insensitive to the precise value assumed for 6 lh . On the other 
hand, at later times (t > 10 days) the radiated luminosity exceeds 
the maximum power due to recombination alone (L > L rec )- This is 
because the thermal energy advected from the base of the wind be- 
comes more important once the mass loss rate (and hence opacity) 
of the wind decreases sufficiently. This transition is evident also in 
the late-time rise of the temperature at the advection radius T adv . 

Figure [8] shown an example of the wind structure with optical 
depth for the M p = lMj case from Figure|9] corresponding to a time 
t = fvisc ~ 1-5 days after disruption when the radiated luminosity 
is L x 1.2 x 10 37 erg s _1 . Note that the ionized fraction X f ss 0.7 
is still relatively high at the trapping radius, but that recombination 
is essentially complete by the time matter reaches the photosphere 
(Xj x 0). Since energy released exterior to the trapping radius es- 
capes without adiabatic losses, this implies that a large fraction of 
the total recombination luminosity L rec (eq. 1151 ) is radiated. 

We conclude that the early super-Eddington accretion phase 
of tidal-disruption planet-star mergers is accompanied by a bright 
optical transient with a peak duration ? visc ~ day-week, which is 
followed by a general dimming and blueward spectral shift lasting 
for a time ?Edd ~ months-year. As we now describe, this dimming is 



Not to be confused with the trapping radius in the disk (eq. 1211 ). 



10 See http : //cdsweb . u- strasbg . fr/topbase/OpacityTables . html 
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Figure 10. Spectra of thermal emission (multi-color blackbody) from 
the sub-Eddington accretion phase (t > t e ^; §5.2\ following the tidal- 
disruption of aplanet, shown at several times r/f e dd = 1 (brown),3 (blue),\0 
(red), and 30 (green), where f c dd =80 days is the Eddington timescale 
(eq. 1241 ). The calculation shown assumes the mass of the disrupted planet 
Mp = lMj and disk viscosity a = 0.1. 

only temporary, since emission re-brightens once radiation from the 
disk itself becomes visible following the super-Eddington phase. 



5.2 Accretion Disk Emission (Sub-Eddington Phase; t > f e( id) 

Powerful outflows from the disk cease once the accretion 
rate decreases below the Eddington rate on a timescale t > 
fedd ~months-year. After this point emission from the accretion 
disk itself becomes visibleFlwith abolometric luminosity L ~ 10 38 
ergs s _1 (M/M c dd.*) that is comparable or greater than the peak lu- 
minosity during the super-Eddington phase (and much greater than 
the luminosity at times / < f cd d; see Fig. [5). The sub-Eddington 
transition thus signals the onset of much brighter and harder emis- 
sion. 

Since the disk is geometrically thin when M < M c dd, its spec- 
trum can be reasonably well described as that of a multi-colored 
blackbody, which extends from the stellar surface r ~ if* to the 
outer edge of the disk r = i?d(0 (eq. 1231 ). The effective tempera- 
ture at any radius in the disk is given by equating the radiated flux 
(TsbTt. with the gravitational power released a (3/87r)GM*Mdr~ 3 . 
The peak frequency of emission resulting from radius r is thus 
given by 

/ M \ 1/4 / T 3/4 
H Mk ,3^ ff ,17ev(-lJ (_) , (29) 

such that we expect the disk emission to range from ultraviolet 
(UV) to optical-infrared wavelengths. 

Figure [10] shows our calculation of the spectrum of thermal 
disk emission for the case of a disrupted planet with mass M p = 
lMj at several times after the Eddington time f c dd = 1, 3, 10, 30 x 
'edd- We again assume a value a = 0.1 for the disk viscosity, which 



In principle ejecta from the super-Eddington phase could still be opaque, 
even long after the outflow subsides. However, a calculation similar to that 
presented in j|4.4l shows that photons from the inner disk will fully ionize 
the outflow, and that only a short period of free expansion is required for 
the ejecta to become optically thin to Thomson scattering. 
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Figure 11. Light curves of thermal emission from the sub-Eddington ac- 
cretion phase at various observed frequencies, including FUV (brown), 
WJV(blue), \J(green), Y(orange), R(purple), and K (red) bands. Top and 
bottom panels show cases corresponding to a disrupted planet of mass 
Mp = lMj and lOMj, respectively. Shown for comparison at times r < ? c dd 
is the visual emission from disk outflows during the preceding super- 
Eddington accretion phase ( jj5.ll see Fig. |5J. The visual light curve dis- 
continuity before and after the Eddington transition indicates that an inter- 
mediate transitional phase (not shown) will smoothly connect these phases. 

results in t eiA ss 80 days (eq. 1241 ). Figure[l0]confirms the expecta- 
tion from equation d29t that the spectrum initially peaks in the far 
UV at t > fedd, but that v pea k(r = Re) K M^ 4 oc r 1/3 decreases 
slowly thereafter (eq. El l29l). 

Figure [TTJ shows the light curves in various wavebands, rang- 
ing from far UV (FUV; v = 2 x 10 15 Hz) and near UV (NUV; 
y = 1.4 x 10 15 Hz) to K-band (v = 1.4 x 10 15 Hz) in the near in- 
frared, shown now for the cases of a planet with mass M p = lMj 
and M p = 10M S in the top and bottom panels, respectively. Observe 
that at optical frequencies, the peak luminosity is vL v > 10 36 erg 
s~' on a timescale > f c dd ~month-year, before fading as vL v oc r 1 
thereafter. This is somewhat greater than the final luminosity from 
the super-Eddington phase, which we have added to the visual light 
curves at times / < f c dd using our results from Figure [9] Compar- 
ing the predicted light curves at t < t eM and t > t M suggests that 
the sub-Eddington transition is probably accompanied by an op- 
tical re-brightening (although the transition could be more subtle 
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depending on the parameters). At UV wavelengths the transition 
is clearly more dramatic, with the flux increasing from effectively 
zero at t <sc f edd to vL v ~ 10 37 erg s -1 at t > t eM . 

Finally, note that our model neglects the ~ 1/2 of the total 
gravitational energy released from the boundary layer connecting 
the disk and star. Radiation from this layer will enhance the UV 
flux by a factor ~ 2 above our predictions in Figures fTOlandfTTI 



6 DISCUSSION AND CONCLUSIONS 
6.1 Summary of Merger Transients 

Table [2] summarizes the predicted electromagnetic transients from 
planet-star mergers described in this paper. 

In the case of direct-impact mergers (ij4} the predicted tran- 
sient is summarized by Figure [7] and quantified in Figure [6] Ini- 
tially, emission is visible directly from the heated stellar surface 
and peaks at EUV/soft X-ray wavelengths. At the earliest times (> 
months-year prior to merger) the signal is periodic on a timescale 
of hours, with a duty cycle ~ 10 per cent and a characteristic lu- 
minosity ~ 10 35 (ex/0.1) erg s (with the exact light curve shape 
depending on the observing angle and geometry of the emitting 
region). Periodicity ends at somewhat later times, but the orbit- 
averaged average luminosity continues to increase exponentially 
with time, before peaking at L EUV / X ~ 10 36 (£x/0. 1) erg s _I on a 
timescale ~week-months prior to the merger. 

EUV/X-ray emission fades once outflows powered by the 
merger becomes sufficiently opaque to block radiation from the 
stellar surface (Fig. [5}. The final stages of the merger are instead 
characterized by optical emission from the outflow itself, which 
brightens from a luminosity L opl ~ 10 35 erg s _1 on a timescale 
of a week-months prior to merger, to a peak luminosity L opt ~ 
10 37 - 10 38 erg s -1 on a timescale of days after merger (eqs. 1171 , 
1161 ). The predicted optical spectrum is that of a dense outflow 
with abundances characteristic of the stellar surface (metallicity 
Z ~ Z Q ). The velocity width is characteristic of the stellar escape 
speed ~ hundreds km s~' and should remain relatively constant in 
time until after the optical peak (after which the line widths may 
decrease as the photosphere recedes to greater depths, and hence 
lower velocities, in the homologously expanding ejecta). 

For tidal-disruption merger events (§5\ Figs. [9l ITOl 1 1 11 1 emis- 
sion initially peaks at optical wavelengths and originates from the 
super-Eddington outflow. This emission brightens on a timescale ~ 
days- week, before peaking at a luminosity L opt ~ 10 37 - 10 38 erg s 
and decaying as vL v cc f" 1 thereafter. At UV wavelengths the emis- 
sion is initially dim, but then brightens considerably to vL Y ~ 10 37 
erg s~' on a timescale / ~ / cdd ~ months- few years, once the hot 
inner edge of the accretion disk becomes directly visible (Fig. II U. 
Spectral line shapes should evolve from P-Cygni profiles charac- 
teristic of a dense outflow at times / <k f cdd , to emission lines char- 
acteristic of the hotter interior at t > f edd . The outflowing material 
likely has a super-solar metallicity (Z ~ 5Z Q ), characteristic of the 
disrupted planet. If the velocity of the outflow is proportional to the 
escape speed from the launching point of the wind, then one pre- 
dicts that the line velocities will initially decrease as the accretion 
disk viscously expands, but then will increase again as the photon- 
trapping radius recedes towards the the stellar surface. Once the 
disk itself becomes visible at times / » / cdd , the emission lines may 
develop double-homed profiles characteristic of the Keplerian disk. 



6.2 Detection Prospects and Strategy 

A planet-star merger in our own galaxy, though intrinsically very 
bright, would be challenging to detect due to extinction through the 
Galactic plane. A more promising detection strategy is via a sys- 
tematic search of nearby galaxies at optical, UV, and X-ray wave- 
lengths. Since to first order the rate of planet-star mergers is pro- 
portional to the total number of stars, the most promising targets for 
such a search are the most massive nearby galaxies, such as M3 1 
(Andromeda). Indeed, since M31 and the Milky Way share a sim- 
ilar mass, one should expect a comparable merger rate, — 0.1 — 1 
yr 1 (©• A systematic survey of M31 over a timescale of years is 
thus one of the most promising ways to detect planet-star mergers. 

One important question is why planetary mergers have not 
yet been detected by existin g surveys. The XMM-Newton survey 
of M31 dPietsch et al1l2005t) identified 18 super-soft X-ray sources 
with luminosities L x > 10 35 erg s" 1 , none of which were obvi- 
ously planet-star mergers. However, since this survey covered a; 1 .2 
square degrees (representing only a fraction of the stellar light of 
M31, which spans several degrees on the sky) and given the low 
predicted rate < 0.1-1 yr~' galaxy -1 of direct-impact mergers, 
it is not particularly surprising that no candidate mergers were de- 
tected. 

A perhaps more promising instrument to search for direct- 
impact mergers is the upcoming satellite mission eROSITA 
(Predehl e t alj|2010h . which has a large field of view ~ 1 degree 
(allowing it to cover M3 1 in just a few pointings) and is sched- 
uled to launch around 2013. A one hour integrated exposure by 
eROSITA would reach a depth at soft X-rays ~ 0.5 - 2 keV of 
~ 10~ 14 erg cm~ 2 s" 1 , corresponding to a luminosity L x ~ 10 36 
erg s _I at the distance of M31, comparable to the predicted peak 
brightness of direct-impact mergers given reasonably optimistic as- 
sumptions (e rad > 0.1). Although the all-sky survey performed by 
eROSITA will not reach this depth for « 4 years (much longer than 
the predicted peak duration of the X-ray emission ~week-months 
from direct-impact mergers), a series of pointed observations over a 
shorter period could reach a substantially greater depth. Candidate 
events could then be followed up with more sensitive observations, 
in order to track the light curve evolution or to search for periodic- 
ity. 

Optical searches also provide a promising means for detect- 
ing planet-star mergers, in particular given the recent advent of 
sensiti ve synoptic surv eys such as the Palo mar Transient F actory 
(PTF: lLawetai] l2009) and Pan-STARRs iKaiser et alj|2002h . Ad- 
vantages of an optical search include the fact that bright emission 
is predicted to accompany both direct-impact and tidal-disruption 
merger events, and that the predicted emission is more robust than 
in X-rays (the latter which depends sensitively on how efficiently 
radiation escapes the stellar surface). The search for "Fast Tran- 
sients in Nearest Galaxies" (P60-FasTING) survey performed with 
PTF dKasliwal et alj[201lh already surveys nearby galaxies with a 
survey cadence of ~ 1 day to magnitude limit g < 21 , correspond- 
ing to a luminosity depth L op , > 6 x 10 36 erg s _I at the distance 
of M31, in principle sufficient to detect the merger of planets with 
masses M p > lMj. 

A major challenge of optical searches is the ability to dis- 
tinguish planet-star mergers from unrelated astrophysical events 
which otherwise appear similar. In particular, optical transients 
from both direct-impact and tidal-disruption mergers may in some 
ways resemble classical novae, which share similar durations and 
peak luminosities (e.g. IShafter et al.l 1201 ll : iKasliwal et alj 1201 ll) 
and are a factor > 50 times more common than even our opti- 
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mistic estimate of the planet-star merger rate. One important dif- 
ference is the characteristic outflow velocity, which is predicted to 
be ~ 100 - 400 km s _1 in the case of planet-star mergers, but is 
typically mu ch higher > 10 3 km s in novae. Of the ~ 40 novae 
compiled by IShafter et alj J201 ll) with measured Ho- line widths, 
only one (M31N 2007-1 lg) had a FWHM velocity < 500 km s" 1 
(see their Fig. 16). Observations of this event one year later, how- 
ever, showed that it had re-brightened to approximately its original 
luminosi ty, indicating that it was probably a long-period Mira-like 
variable dShafter e t al. 2008, ATel 185 1) instead of a nova or planet 
merger event. 

Another distinguishing feature of novae is that they result from 
the ejection of a discrete shell of ejecta (powered by a thermonu- 
clear runaway), followed by a relatively steady wind from the sur- 
face of the hot white dwarf. For direct-impact mergers (^4), by 
contrast, we predict the opposite evolution- namely, optical emis- 
sion that originates from a quasi steady-state wind at early times, 
before later transiting into a dynamical ejection just prior to the 
merger. Tidal-disruption events (^5} in some ways represent a 
closer analog to novae, in that the event evolves fastest initially, 
after which the 'hot' central object (the inner accretion disk in the 
case of planet-star mergers, versus the hot white dwarf in novae) 
becomes increasingly visible with time. The predicted light curve 
in this case (rise time ~ days-week, followed by a subsequent de- 
cay vL v oc r 1 ) also bears a qual i tative resemblance to that of some 
novae (e.g. lKasliwal et alj|201 lklShafter et alj201 ll) . 

Ultimately, the most promising strategy to detect and identify 
planet-star mergers is to combine optical, UV, and X-ray obser- 
vations. For instance, one way to distinguish novae from direct- 
impact planet mergers is to search for spatially coincident X-ray 
emission months-year prior to optical maximum. Unfortunately, 
most searc hes to date have focused on X-ray e mission following 
the nova jPietsch et al .11 20051 : iHenze et alj|20ld ). with only a few 
searc hes for 'quiescent' emission prior to eruption (e.g. lOrio et al.1 
1200 ll) . Likewise, a promising way to distinguish novae from planet 
tidal-disruption events is to search for spatially coincident bright 
UV emission months-years following optical maximum. Multi- 
wavelength diagnostics, aided by systematic surveys of nearby 
galaxies, are clearly key to successfully identifying these relatively 
rare events. 



6.3 Conclusions and Future Work 

Basic considerations suggest that planets merge with their host stars 
at a relatively regular rate (^2j Fig.[TJi. Detecting such events di- 
rectly, beyond confirming the basic physical picture outlined here, 
could provide key information on the properties of the merging bi- 
nary, as well as inform the rate of tidal dissipation (i.e. the value of 

and/or the rate at which planet s currently flow into their inner 
stellar system dSocrates etal. 2011). 

In §[3] we identified three qualitatively different merger out- 
comes ('direct-impact', 'tidal-disruption' and 'stable mass trans- 
fer'), which, from the current sample of explanets, we estimate that 
an appreciable fraction of systems will experience each (Figs.[2j 
Table[T]l. Our estimate of the relative fraction of merger outcomes 
assumed that the current masses and radii of extrasolar planets 
and their host stars represent their values at the point of merger, 
an assumption of questionable validity. Although in ^3] we de- 
scribed the effects of late stellar evolution on the frequency and 
outcome of planet-star mergers, the properties of the planet itself 
may change even when the star is still on the main sequence. As 
a planet migrates inwards towards merger, it experiences increased 



interior heating (e.g. tidal dissipation, ohmic heating, atmospheric 
irradiation) which may inflate its radius over the nominal value 
given its age (e.g. iBodenheimer et al1l200ll ; iBurrows et al. I l2008l : 
llbgui & Burrows 2009k iBatygin et al]|201lh . Close-in planets may 
furthermore lose mass to outflows powered by stellar irradiation, 
although current models suggest that ho t Jupiters are unlikely to 
completely e vaporate by this process (e.g. lMurrav-Clav et alj2009l 
lAdamsll201 ]]; thoug h less massive planets may be more susceptible 
to evaporation, e.g. iRappaport et al.ll2012h . Given these uncertain- 
ties, it is difficult to determine with confidence which merger out- 
come is the most frequent. 

In this paper we have focused primarily on the merger of plan- 
ets in quasi-circular orbits, since most known hot Jupiters should 
merge with final eccentricities e < 10~ 3 (Appendixf/jt. Such a small 
eccentricity would not qualitatively change our conclusions regard- 
ing the dynamics of inspiral described in 94. II since the difference 
between pericenter and apocenter is smaller than the stellar scale 
height H > lO^R*. Nevertheless, depending on the rate at which 
planets are deposited directly into the inner stellar system by planet 
scattering or Kozai oscillations (or to what extent gravitational in- 
teractions with other planets pump eccentricity), some planets may 
enter with low impact parameters on highly eccentric orbits. 

Depending on the impact parameter of the trajectory and the 
mass ratio of the system, an eccentric planet-star merger can result 
in either a direct head-on co llision; complete eject ion of the planet; 
or a tidal-disruption event (iGuillocho n et al.ll201ll) . In the case of 
a direct head-on collision, the transient may resemble a shortened 
version (or final stages) of a direct-impact merger (i}4j, possibly 
producing a prompt flash of high-energy emiss ion from the stel- 
lar surface (see also Zhang & Sigurdsson 2003) or thermal emis- 
sion from an impact-driven outflow ( 34.4.2^ In the tidal-disruption 
case, the emission would qualitatively resemble the quasi-circular 
merger case that we have already discussed (jj3), except that the 
late-time accretion rate (and hence the resulting optical-EUV light 
curve) wil l decrease d ifferently with time M oc r 5/3 due to fall-back 



accretion < lReeslll988l) . than the evolution M 



oc t 



-4/3 



predicted due 



to the viscous evolution of the disk (eq. 1221 ). 

In addition to classical novae, transient emission from planet- 
star mergers may in some ways resemble that from the merger 
between two binary stars, events which may occur at a simi- 
lar rate. The optical transients M85 OT2006-1 an d V838 Mon 
have b oth been argued to result from s t ellar merge rs (Muna rTet al.l 
120021 ; ISoker & Tvlendal 120031 l2006t iKulkarni et al.l 120071) . al- 
thou gh other explanat i ons have been proposed for these events 
(e.g. iRetter & Maroml 120031 iPastorello etall 120071) . The recent 
OGLE transient VI 309 Scorpii repre sents an essentially indis- 
putable case of a binary stellar merger (lTvlendaetal.ll201ll) . 

One might think that a merger between two stars would be 
substantially brighter than a planet-star merger due to the consid- 
erably greater gravitational energy released during such an event. 
However, this expectation may not necessarily be realized. One 
"advantage" enjoyed by low mass-ratio mergers is that insufficient 
energy is released during the earlies stages of the merger to ap- 
preciably expand the outer stellar envelope. In this case the heated 
atmosphere of the star may have sufficient time to cool before the 
next planetary orbit, thus increasing the likelihood that high energy 
emission may escape the viscnity of the stellar surface ( 94.4. It . at 
least during the earliest stages of the merger. 

A perhaps more important advantage of planetary versus stel- 
lar mergers is that, in binaries with mass ratios q > 0. 1 the sec- 
ondary is sufficiently massive to maintain synchronous rotation 
of the stellar envelope of the primary until the point of merger 
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Table 2. Key Properties of Transients from Planet-Star Mergers 



Event Type 


Photon Energy 


Peak Luminosity 
(erg s _1 ) 


Peak Duration 


Ejecta Metallicity 


Characteristic Velocity 


Direct-Impact (fl Figs.[3]|6|[7): 


Stellar Surface Emission (34.4. 1> 
Inspiral-Driven Outflow (34.4.21 


EUV/Soft X-ray 
Optical 


10 36 fe ad /0.1) 
10" - 10 38 (eq. (HI) 


months-year 
days (eq. 1161) 


N/A 

~Z G 


N/A 

~ 100 - 400 km s" 1 (blue shifted) 


Tidal-Disruption (@ Figs.l9l[T0lfni: 


Super- Eddington Disk Wind (35.lt 
Accretion Disk (35.21 


Optical 
Optical(UV) 


10 36 - 10 37 (M P /Mj) 
10 36 (10 37 ) 


day-week 
months -year 


~5Z 
N/A 


~ 100 - 400 km s" 1 (blue shifted) 
~ 100 - 400 km s" 1 (double horned) 


Classical Novae*"' 


Optical 


10 37 - 10 39 


days-months 


~ 3 - 30Z* 1 


> 10 3 km s _1 (blue shifted) 



™ Shown for purposes of comparison; (fc lLivio & TruranH 19941) . 



(Soker & Tvlenda 2006]). This implies that substantially less rela- 
tive kinetic energy is available between the merging stellar surfaces 
to pow er prompt emission than in the case of a less massive sec- 
ondary llSoker & Tvle nda 2006). The quantity of mass ejected, and 
the resulting prompt electromagnetic display, could thus in prin- 
ciple be much less luminous than one would anticipate based on a 
simple 'scaled-up' version of the model presented in 94.4l (in which 
we have implicitly assumed a slowly rotating star). 

Despite these caveats, to the extent that stellar mergers eject 
some quantity of ionized material, much of the same physics de- 
scribed in ^4] and §|5] regarding optical transient emission should 
apply to these events as well. For instance, VI 309 reached a peak 
optical luminosity ~ 10 38 erg s -1 lasting for a timescale ~ 1 month, 
such that the total radiated energy was ~ 3 x 10 44 erg. If this lu- 
minosity were powered primarily by the recombination of hydro- 
gen in an outflow (as we find characterizes our outflow models in 
94.4.21 and 95.1b . then the requisite mass of ejecta is ~ 0.01M o , a 
quantity in agreement w ith independent estimates of similar events 
dTvlenda & Soke3l2006l) . 

Much additional work is required to better understand the ob- 
servable signatures of planet-star mergers. In the case of direct- 
impact mergers, a more detailed model of the hot post-planetary 
layer ( 94.2b is necessary to make more quantitative predictions, the 
study of which will likely require a detailed radiation hydrodynam- 
ical calculation. Our basic estimates hint that a sizable fraction of 
the thermal energy generated near the surface may escape as high 
energy radiation; however, a better understanding of the detailed 
structure of this region (including the effects of sustained heating 
in the vicinity of the planet-star interface) is necessary to assess 
the radiative efficiency e ra[ j, geometry, and beaming fraction of the 
resulting emission (Fig. [3}- Another major uncertainty is the ex- 
tent to which the heated stellar surface is able to cool, or to what 
extent the star remains 'puffed up', by the time the planet has com- 
pletely orbited the star. Our calculations in 94.11 assumed that the 
planet interacts with an essentially unperturbed stellar atmosphere 
after each passage, clearly an oversimplification, especially during 
the final dynamical stages of the merger. Depending on the detailed 
structure of the perturbed atmosphere, this could accelerate (or de- 
celerate) over the inspiral evolution calculated in 94.11 

Given the sensitive, wide-field capabilities of surveys at X-ray 
and optical wavelengths, it seems promising that planet-star merg- 
ers should be detected within the next decade. 
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APPENDIX A: STABLE VS. UNSTABLE MASS 
TRANSFER 

Here we quantify the fate of a planet that fills its Roche-lobe outside 
the stellar surface, i.e., for which a, > R+ + X, (eq. (3)), where 
a, (eq. (2)) is the Roche radius and X, is the distance of the LI 
Lagrange point from the center of mass as described in ^3] Our goal 
is to explore whether infinitesimal mass loss from the planet dM p < 
through the Lagrange LI point results in the planet expanding 
faster or slower than the size of Roche-lobe, in order to determine 
whether mass transfer is stable or uns tableful 

Equating the specific angular momentum j m = 
(GMJaj'^fa, - X,) of matter leaving the planet with that 
of a Keplerian orbit about the star = (GM+R c ) l/2 we obtain the 
characteristic circularization radius of the material 

R c =c h (l-X,la,f. (Al) 

If R z < R t then the accretion stream directly impacts the stel- 
lar surface and no disk forms, such that all of the mass and angular 
momentum carried by dM p is transferred directly to the star. In the 
much more common case that a disk forms (i.e. R c > R+) it is 
still reasonable to assume that most of the transferred mass will 



12 For the extreme mass ratio systems under consideration, the distance 
from th e center of the planet to the L2 Lag range point is similar to that to 
LI (e.g. lGuetalHc XB; Nayakshin & Lodato 201 fj). Some mass will there- 
fore overflow from L2 as well. Mass loss from LI is, however, still likely to 
dominate since the star is strongly irradiated from the front, resulting in a 
larger atmospheric scaleheight near LI. Since mass loss from L2 results in 
angular momentum loss from the system, its effect is to increase the likeli- 
hood of unstable mass transfer and dynamical tidal-disruption (see below). 
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ultimately be accreted the star(3 however, the fate of the angular 
momentum in this case is less clear. 

Assuming that the disk is Keplerian, its specific angular mo- 
mentum near the stellar surface is j im = (GM^i?*) 1 ' 2 . Thus only a 
fraction / = j lm / j m = (R^/Rc)' 12 of the angular momentum of the 
accreting material is transferred onto the star. The rest remains in 
the disk, where it will eventually be returned to the orbit via tides 
exerted on the planet. 

Mass loss reduces the angular momentum of the planet by 
an amount dj p = dM p (GM+a) i/2 + M P (GM* /a) l/2 da/2. Equating 



this to the angular momentum gained by the star dj p 
dM p (GM+R+) l/2 results in the following expression 



dJ„ h 



— = -2— h - (— V /2 1 

dM„ Mj \ a I ^ 



The volume of the Roche-lobe scales as V RL 
dV RL 



oc M p a 3 or 



M„ aAMJ Mr, [ \ ^ a > >\ 



(A2) 



(A3) 



1 p lrl p 

For a = a, (eq. |2)) the bracketed quantity in the last expression 
takes the value [3 V2(pp/p*)' /6 - 5], Equation dA3l > shows that if 
overflow occurs at large orbital separations (a, » R+ or p p /p* << 
1), then the volume of the Roche-lobe increases rapidly with mass 
loss, while for a, ~ R+ its size evolves more gradually. 

The change in the radius of a planet that undergoes mass 
loss depends on its entropy and ma s s, i.e. R P (M P ,S). The plan- 
etary models of Spie gel & Burrows! d2012h suggest that the rel- 
atively large observed radii of close-in planets indicate high en- 
tropies S ~ %fi/baryon. If one assumes that the planet loses mass 
at consta nt entropy, then the analytic fits to the mass-radius relation 
given by Spiegel & Burrows (2012) give the following expression: 



dV„ 



C- 



(A4) 



where C = 31nl0[(p m + puS + p 2l S 2 )(M p /M,) + 2(p 02 + 
Pi 2 S)(M D / 'Mj) 2 ]; S is in units of fc B /bar yon; and the parameters 
Pij are given in Spiegel & Burrows (2012). Since typical values are 
C ~ -1 for M p ~ lMj and C ~ -0.3 for M p ~ 10Mj, we conclude 
that most of the observed planets expand upon mass loss. 

Since the masses of the planets in our sample are measured, 
we can directly estimate the value of j^- for individual systems 
that will undergo Roche-lobe in order to predict their ultimate fate. 
In particular, if \-^-\ > \^^\ then the planet will expand faster 
than the Roche-lobe with mass loss, indicating that mass transfer 
is unstable (and vice versus). Figure [2] summarizes our results ap- 
plying this analysis to infer the stability of mass transfer for the 
entire sample of transiting planets (assuming that the properties of 
the planet and star today reflect those at the time of merger). 



APPENDIX B: JUSTIFICATION OF CIRCULAR ORBITS 
FOR DIRECT-IMPACT MERGERS 

Here we analyze the tidal evolution of star-planet systems in order 
to derive the orbital decay timescale and to estimate the eccentricity 
of systems just prior to merger or tidal-disruption. We focus on 
a compact planet-star systems and ignore the effects of potential 
interactions with other planets (e.g. eccentricity 'pumping'). 

13 Although this n eed not always be the case, especially in the case o f ec- 
centric disruptions l lGuillochon et al.l201 ltlkavakshin & Lodatol201 3) . 



Assuming that the planet resides in a synchronous orbit 
about a relatively slowly rotating star, the evolution of the or- 
bital parameters of the system e, a are giv en by the expres- 
sions ( in the l ow ec c entricity e < 0-2 limit; see Goldreich & SoteJ 
I l966l: iKaulal 1 19681: IPeale & Casserj 1 19781 : Ijackson et alj |200S| ; 
barnes et alj2009l:llbgui et alfeplOh 



-13/2 



[KiB*/QL + K 2 R 5 JQ'J, 



(Bl) 



1 de 

e At 

1 da = _ a -W2 [2KiR s 2 /Q , + SK 2 R 5 J25Q'J, (B2) 
a at 

where K { = (63 / 4)G 1/2 Ml' 2 /M p , K 2 = (225 /16)G' 12 M p /M l J 2 , and 
Q' p is the modified tidal quality factor of the planet; the rest of 
the quantities are defined in Section ^3] For a discussion of the 
several underlyin g assumptions in deriving these expressions, see 
llbgui et aljj2O10h . 

Tides raised on the planet (terms containing Q' p ) serve to cir- 
cularize the orbit. For sufficiently low Q' p , the orbit circularizes on a 
timescale shorter than that of orbital decay. The fact that the shorter 
period planets with P < 3 days are in near circular orbits (e < 0. 1, 
in contrast to wider systems that can be very eccentric) indicates 
that circulari zarion tides operate effectively in the tight star-planet 
systems fe.g. Lfackson et alj2 008). If this is the case, then both tidal- 
disruptions and mergers take place when the orbit is nearly circular. 

An upper limit on the final eccentricity of the orbit is derived 
by setting Q' p — > <x>, in which limit a and e decay on the same 
timescale. The ratio of expressions B1/B2 in this limit leads to 
de/e = (25/8)da/a. Integrating the last expression starting from 
the current values a iu , e m one finds e f = e m (a f /a iu ) 2S/>s . The most 
compact systems with a m ~ (6 — S)R Q have observed e m < 0. 1. Set- 
ting fl f = R Q (merger case) gives e f < 3 x 10~ 4 while for a { = 2R Q 
the final eccentricity of the orbit is e { < 3 x 10~ 3 . These small val- 
ues of eccentricity at the point of tidal-disruption or merger are not 
expected to appreciably change the analysis presented in §0 which 
implicitly assumed perfectly circular orbits. This is because the dif- 
ference between pericenter and apocenter ~ 2tf < 10~ 3 is smaller 
than the stellar scale height H > 10~ 3 R*. 

Exression B2 also provides an estimate for the orbit tidal de- 
cay timescale (again setting Q' p — » oo) as given in equation (TJ 
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